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ABSTRACT 

Nineteen of the ~40 RR Lyrae stars in the Kepler field have been identified as can- 
didate non-Blazhko (or unmodulated) stars. In this paper we present the results of 
Fourier decomposition of the time-series photometry of these stars acquired during the 
first 417 days of operation (Q0-Q5) of the if epZer telescope. Fourier parameters based 
on ^18400 long-cadence observations per star (and ^150000 short-cadence observa- 
tions for FN Lyr and for AW Dra) are derived. None of the stars shows the recently 
discovered 'period-doubling' effect seen in Blazhko variables; however, KIC 7021124 
has been found to pulsate simultaneously in the fundamental and second overtone 
modes with a period ratio P2/P0 ~ 0.59305 and is similar to the double-mode star 
V350 Lyr. Period change rates are derived from O— C diagrams spanning, in some 
cases, over 100 years; these are compared with high-precision periods derived from the 
Kepler data alone. Extant Fourier correlations by Kovacs, Jurcsik et al. (with minor 
transformations from the V to the Kp passband) have been used to derive underlying 
physical characteristics for all the stars. This procedure seems to be validated through 
comparisons of the Kepler variables with galactic and LMC RR Lyrae stars. The most 
metal-poor star in the sample is NR Lyr, with [Fe/H]= —2.3 dex; and the four most 
metal-rich stars have [Fe/H] ranging from —0.6 to -1-0.1 dex. Pulsational luminosities 
and masses are found to be systematically smaller than L and M values derived from 
stellar evolution models, and are favoured over the evolutionary values when periods 
are computed with the Warsaw linear hydrodynamics code. Finally, the Fourier pa- 
rameters are compared with theoretical values derived using the Warsaw non-linear 
convective pulsation code. 
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1 INTRODUCTION 

The Kepler Mission is designed to detect {via transits) 
Earth-like planets around solar-type stars (Koch et al. 
2010). To achieve this goal the 1.4-ni Kepler telescope has 
been monitoring almost continuously the light variations of 
over 150 000 stars in a field located at (RA,DEC)j2ooo = 
(19:22:40, -^44:30) or (Z,fo) = (76.3, -hl3.5). The mandate of 
the Kepler Asteroseismic Science Consortium (KASC) is to 
use these data to better understand the astrophysics of the 
stars seen in the Kepler field, with Working Group 13 being 
responsible for the study of the RR Lyr variable stars. 

Kolenberg et al. (2010) discussed the initial selection of 
28 target RR Lyr stars in the Kepler field and presented pre- 
liminary results based on the earliest Kepler observations for 
RR Lyr and V783 Cyg, both of which exhibit the Blazhko 
(1907) effect, and for the 'non-Blazhko' star NR Lyr. An im- 
portant discovery was the phenomenon of 'period doubling', 
a cycle-to-cycle amplitude variation that occurs at certain 
Blazhko phases. The effect has been explained by Szabo et 
al. (2010) as a 9:2 resonance between the fundamental mode 
and the 9th-order radial overtone. More recently KoUath, 
Molnar & Szabo (2011) showed that the resonance desta- 
bilises the fundamental mode limit cycle leading to period 
doubling. Among RR Lyr stars period doubling has been 
observed to date only in the Blazhko variables. Benko et 
al. (2010, hereafter BIO) studied a sample of 29 modulated 
and unmodulated RR Lyr stars using Kepler Q0-Q2 long- 
cadence data, and identified 14 of the 29 stars (48%) as ex- 
hibiting the Blazhko effect. Most recently, a more detailed 
investigation of RR Lyrae itself has been made by Kolenberg 
et al. (2011), and Smolec et al. (2011) have used non-linear 
hydrodynamic models to model Blazhko RR Lyr stars. 

In the present paper we look more closely at the unmod- 
ulated ab-type RR Lyr stars - that is, those stars pulsating 
in the fundamental mode and believed to exhibit no ampli- 
tude or phase modulations and to have the most stable light 
curves. According to the notation established by Alcock et 
al. (2000) our sample stars are of type RRO. The analy- 
sis is based primarily on the Kepler long-cadence (30-min) 
photometry acquired in quarters Q0-Q5 (over 18000 mea- 
surements per star made over 417 d), and the short-cadence 
(1-min) photometry for FN Lyr and AW Dra acquired in 
QO (over 14000 high-precision brightness measurements per 
star made over 9.7 d) and Q5 (~135,000 observations per 
star made over ~90 daysJJ- 

We first discuss the sample selection (§2.1) and proce- 
dure for transforming the raw fluxes delivered from the Ke- 
pler Science Office to magnitudes on the Kp system (§2.2). 
Before analyzing the Kepler data a search (§3.1) was con- 
ducted of available historical data for the stars [e.g., pho- 
tometry, periods, light curves, and times of maximum light). 
Accurate periods were derived using all the available Kp 
photometry (§3.2). Then, from the plotted light curves, to- 
tal amplitudes and risetimes were measured (§3.3). Period 

^ Public data are available from the internet website 
: //archive. stsci.edu/kepler/l. It is the policy of the Kepler 
Mission to make public all the data mentioned in published re- 
search papers. Thus the Q0-Q5 long and short cadence data anal- 
ysed here will become publicly available upon publication of the 
present paper. 



change rates {AP/At) and revised periods were calculated 
from O— C diagrams constructed with the historical data 
combined with the Kepler data (§3.4). Fourier analysis of 
light curves (§4.1) was carried out for both the Kp data 
(§4.2) and for the extant V data (§4.4). The long time 
base and intensive sampling of the Kepler data provided 
an opportunity to examine the light curves for cycle-to- 
cycle variations. To test for such variations the data were 
divided into single-cycle blocks, and for each block a sepa- 
rate Fourier analysis was conducted (§4.3). For the Kepler 
stars observed correlations among the Fourier parameters 
(all cycles combined) were examined (§5) for both the Kp 
and V passbands (§5.1). The resulting V-Kp offsets (§5.2) 
allowed comparisons to be made between the Kepler stars 
and RR Lyr stars in galactic and Magellanic Cloud globu- 
lar clusters (§5.3) and field stars in the inner regions of the 
Large Magellanic Cloud (§5.4). Physical characteristics were 
also derived from suitably-modified extant V-band correla- 
tions (§6). These include iron (or metal) abundance [Fe/H], 
dereddened colour (B — V)o, effective temperature TcS, dis- 
tance d, absolute magnitude Mv , pulsational and evolution- 
ary luminosities, L(puls) and L(evol), pulsational and evo- 
lutionary masses, 7Vf(puls) and M{evol), and location in 
the instability strip. Finally, new non-linear hydrodynami- 
cal models (Warsaw code) are presented and compared with 
the observational data (§7), and our results are summarized 
(§8). In Appendix A the newly discovered double- mode star 
KIC 7021124 is discussed, and in Appendix B the cyclic be- 
haviour of nine non-Blazhko stars in the {V — I, l/)-diagram 
(based on ASAS-North data) is discussed. 



2 DATA PROCESSING 

2.1 Selection of non-Blazhko stars 

Approximately 40 RR Lyr stars are presently known in the 
Kepler field, up from the 29 stars discussed by Kolenberg et 
al. (2010) and BIO. Four are multiperiodic c-type RR Lyr 
stars and are the subject of another paper (Moskalik et al. 
2011, in preparation), and ~20 are Blazhko stars, many of 
which have quite complex modulated light curves (see BIO). 
We are interested here in the unmodulated non-Blazhko ab- 
type RR Lyr stars with the most stationary light curves. 
When the Fourier parameters and physical characteristics 
of these stars are known this information will provide a use- 
ful reference for future comparisons with the more complex 
Blazhko stars (see Jurcsik et al. 2009). 

Table 1 of BIO gives coordinates for 16 of the 19 stars 
considered here, and pulsation periods to five decimal places. 
Three stars are not in the BIO: FN Lyr (KIC 6936115) and 
AW Dra (KIC 11802860), both of which were observed at 
short-cadence (SC) in QO with long-cadence Ql data and 
short-cadence Q5 data released later; and KIC 7021124, 
which was discovered subsequently. The RA and DEC 
coordinates (J2000) of these three stars are as follows: 
FN Lyr (19:10:22.25, +42:27:31.6), AW Dra (19:00:48.00, 
-f50:05:31.3) and KIC 7021124 (19:10:26.69, -f 42:33:37.0). 

Three of the stars included in the present study are spe- 
cial cases: (1) V349 Lyr was identified by BIO as a Blazhko 
star with a small amplitude modulation and Blazhko-period 
longer than 127 days. It is one of the faintest RR Lyr stars 
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Table 1. Basic data for the 19 Kepler non-Blazhko ab-type RR Lyr stars 



KIC Star {Kp) Period [d] Freq. [d"^] to RT{Kp) Kp-range AtotiKp) 

name [mag] (this paper) (this paper) [BJD] [phase] [mag] [ma^g] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



3733346 


NR Lyr 


12, 


.684 


0, 


.6820264(2) 


1.4662189(4) 


54964, 


,7381 


0, 


,144 


12, 


,215- 


■12, 


.982 


0, 


.767 


3866709 


V715 Cyg 


16, 


.265 


0, 


.47070494(4) 


2.1244731(2) 


54964, 


,6037 


0, 


,136 


15, 


,612- 


■16, 


.600 


0, 


.988 


5299596 


V782 Cyg 


15 


.392 


0, 


.5236377(1) 


1.9097173(4) 


54964. 


,5059 


0, 


,200 


15, 


,108- 


■15, 


.631 


0, 


.523 


6070714 


V784 Cyg 


15 


.370 


0, 


.5340941(1) 


1.8723292(4) 


54964. 


,8067 


0, 


,195 


15, 


,036- 


■15, 


.670 


0, 


.634 


6100702 


KIC 6100702 


13, 


.458 


0, 


.4881457(2) 


2.0485687(8) 


54953, 


,8399 


0, 


,200 


13, 


,140- 


■13, 


.715 


0, 


.575 


6763132 


NQ Lyr 


13, 


.075 


0, 


.5877887(1) 


1.7012916(3) 


54954, 


,0702 


0, 


,144 


12, 


,576- 


■13, 


.387 


0, 


.811 


6936115 


FN Lyr 


12 


.876 


0, 


.527398471(4) 


1.89609954(1) 


54953. 


,2690 


0, 


,118 


12, 


,146- 


■13, 


.227 


1, 


.081 


7021124 


KIC 7021124 


13 


.550 


0, 


.6224925(7) 


1.606445(2) 


54965, 


,6471 


0, 


,139 


13, 


,018- 


■13, 


.849 


0, 


.831 


7030715 


KIC 7030715 


13, 


.452 


0, 


.6836137(2) 


1.4628145(4) 


54953, 


,8434 


0, 


,177 


13, 


,092- 


■13, 


.739 


0, 


.647 


7176080 


V349 Lyr 


17, 


.433 


0, 


.5070740(2) 


1.9720987(8) 


54964, 


,9555 


0, 


,130 


16, 


,780- 


■17, 


.768 


0, 


.988 


7742534 


V368 Lyr 


16 


.002 


0, 


.4564851(1) 


2.1906520(5) 


54964. 


,7828 


0, 


,134 


15, 


,237- 


■16, 


.370 


1, 


.133 


7988343 


V1510 Cyg 


14 


.494 


0, 


.5811436(1) 


1.7207451(3) 


54964. 


,6695 


0, 


,138 


13, 


,868- 


■14, 


.848 


0, 


.980 


8344381 


V346 Lyr 


16, 


.421 


0, 


.5768288(1) 


1.7336166(3) 


54964, 


,9211 


0, 


,141 


15, 


,809- 


■16, 


.773 


0, 


.964 


9508655 


V350 Lyr 


15, 


.696 


0, 


.5942369(1) 


1.6828305(3) 


54964, 


,7795 


0, 


,143 


15, 


,086- 


■16, 


.059 


0, 


.973 


9591503 


V894 Cyg 


13 


.293 


0, 


.5713866(2) 


1.7501285(6) 


54953. 


,5627 


0, 


,144 


12, 


,579- 


■13, 


.684 


1, 


.105 


9947026 


V2470 Cyg 


13 


.300 


0, 


.5485894(1) 


1.8228569(3) 


54953. 


,7808 


0, 


,193 


12, 


,971- 


■13, 


.570 


0, 


.599 


10136240 


V1107 Cyg 


15, 


.648 


0, 


.5657781(1) 


1.7674774(3) 


54964, 


,7532 


0, 


151 


15, 


151- 


■15, 


.969 


0, 


.818 


10789273 


V838 Cyg 


13, 


.770 


0, 


.4802799(1) 


2.0821192(4) 


54964, 


,5731 


0, 


134 


13, 


,040- 


■14, 


.138 


1, 


.098 


11802860 


AW Dra 


13, 


.053 


0, 


.6872160(2) 


1.4551466(6) 


54954. 


,2160 


0, 


165 


12, 


,538- 


■13, 


,430 


0, 


.892 



in the Kepler field, at (Kp) — 17.433 mag. The detrended 
light curve in fig. 1 of BIO shows the amplitude decreas- 
ing slowly throughout the Q1-Q2 epochs. Classification of 
this star depends on the detrending procedure that is used 
to transform the raw fluxes to magnitudes (see next sec- 
tion and Fig. 1). It has been included here as a borderline 
non-Blazhko star. (2) V350 Lyr was discovered by BIO to 
be pulsating simultaneously in the fundamental and second 
overtone modes (i.e., an RR02 star), with a secondary pe- 
riod P2 = 0.35209 d and a period ratio P2/P0 = 0.592. So 
far there has been no evidence for the Blazhko effect. With 
the additional ~270 days provided by the Q3-Q5 observa- 
tions this star's behaviour has been conflrmed and strength- 
ened. One question of interest concerns the stability of the 
/2 frequency - BIO determined that the /2 amplitude was 
varying between 0.00 and 0.05 mag. This result is conflrmed 
by the newer data. Recent work to better understand this 
phenomenon has been presented by KoUath et al. (2011). (3) 
KIC 7021124 was found during the course of the present 
study to be pulsating simultaneously in the fundamental and 
second-overtone modes, and to be similar to V350 Lyr. The 
details of its pulsational behaviour are given in Appendix A. 



Basic data for the 19 unmodulated Kepler RR Lyrae 
stars is given in Table 1. Column 3 contains the mean 
Kp magnitude given in the Kepler Input Catalog (KIC). 
Columns 4-5 contain our best estimate of the pulsation pe- 
riod and frequency derived using the Kp data (see Table 3 
for details). Column 6 contains accurate values for the times 
of maximum fight, and column 7 lists precise risetimes (RTs) 
derived from the ifp-photometry. Columns 8-9 contain the 
Kp magnitudes at which the light curve reaches maximum 
and minimum fight, the difference being the total amplitude 
in the Kp passband. 



2.2 From raw fluxes to Kp magnitudes 

The preprocessing phases of the long-cadence (LC) Kepler 
data have been described by Jenkins et al. (2010a,b) and 
BIO. Since then three more 'seasons' (or ~270 days) of LC 
data (Q3-Q5) have become available for 18 of the 19 stars 
in our sample. Also available now are QO (~10 days) and 
Q5 (~90 days) SC photometry (see Gilliland et al. 2010) for 
FN Lyr and AW Dra. The procedures used here for trans- 
forming the raw fluxes to Kp magnitudes are described next. 

2.2.1 Long-cadence data 

The present analysis of the stars in Table 1 is based on the 
raw fluxes. Sixteen of the stars were observed at LC (30 min) 
in 'seasons' Ql to Q5, and of these, flve were also observed in 
QO. KIC 6100702 was not observed in Q5, and KIC 7021124 
was observed only in Ql. The Q2-Q4 LC observations of 
FN Lyr were acquired as part of the Kepler Guest Observer 
Program and Drs. K. Mighell and S. Howell kindly shared 
their observations with us. The original data consisted of 
raw flux counts as a function of time, with times given as 
barycentric Julian date (BJD). For convenience all the times 
were translated by subtracting 54953 d so that the flrst QO 
observations occur at time t = 0.53 d. 

Figure 1 helps to illustrate the procedures used to con- 
vert long cadence fluxes to Kp magnitudes. The left pan- 
els show the observed raw fluxes (loglO scale) of V894 Cyg 
(top) and V349 Lyr (bottom) plotted against time. Both 
stars were observed every 30 minutes in 'quarters' Ql to Q5 
and V894 Lyr was also observed in QO. One sees zero point 
differences (sensitivity variations) resulting from quarterly 
'rolls' of the Kepler telescope causing the light from a given 
star to fall on different CCD chips and pixels each quar- 
ter, and trends (both linear and non-linear, and increas- 
ing and decreasing) within a given quarter. The scalloped 
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Va94 Cyg (KIC 9691503, <J(:p>= 13.393, LC:Q0-Q5) 
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\ \ \ 
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V349 Lyr (KIC 7176080, <{fp>= 17.433, Ql-QB) 



100 300 300 

BJD - 54953 
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300 300 
BJD - 64963 



Kp 



V349 Lyr (KIC 7176080, Kp=17-433, Q1-Q5) 



100 300 

BJD - 64963 



Figure 1. Raw (left) and processed (right) pliotometry for V894 Cyg (top) and V349 Lyr (bottom). V894 Cyg was observed 18911 
times at long eadence in quarters Q0-Q5 and is of intermediate brightness with raw flux density ~ 1.6 X 10* e~ /cadence. V349 Lyr is 
the faintest RR Lyr star in the ft" epier field and was observed 18387 times at LC:Q1-Q5 with raw flux density ~ 2 X 10^ e^ /cadence. 



pattern seen for V894 Cyg (with amplitude larger at max- 
imum than at minimum light) is due to the 30-min inter- 
val between the LC observations going in and out of phase 
with the light cycle determined by the pulsation period of 
0.57 d. The time ranges for each LC 'season' were as fol- 
lows: BJD 54953.53 to 54963.24, corresponding to time 0.53 
to 10.24 d (QO), BJD 54964.51 to 54997.98, or 11.51 to 
44.98 d (Ql), BJD 55002.94 to 55091.47, or 49.94 to 138.47 d 
(Q2), BJD 5093.52 to 55182.50, or 140.22 to 229.50 d (Q3), 
BJD 55185.38 to 55275.20, or 232.38 to 322.20 d (Q4), and 
BJD 55275.99 to 55370.66, or 322.99 to 417.66 d (Q5). 

The right panels of Fig. 1 show the corresponding de- 
trended and normalized Kp apparent magnitudes. The pro- 
cedure for transforming the raw fluxes into apparent magni- 



tudes was as follows: (1) the data were separated into blocks, 
usually by 'quarter' but various blocks were tried (see the 
left panels of Fig. 1); (2) the raw fluxes were converted to 
apparent magnitudes using 'mag = 33.5 — 2.5xloglO(Fraw)', 
where the magnitude zero point is as yet only approximate 
but gives mean values roughly in accord with the KIC mean 
magnitudes; (3) within a block the magnitudes were de- 
trended by fitting and subtracting either a linear or polyno- 
mial fit; (4) the detrended data for each block were shifted 
to the Kp mean magnitude level given in the KIC; (5) within 
each block the magnitudes were stretched or compressed by 
multiplying by an appropriate scale factor so that the mini- 
mum and maximum magnitudes were in agreement with an 
adopted reference quarter (usually either Q4 or Q5); (6) ob- 
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vious outliers (cr > 5) were removed - before this was done 
the total number of LC observations per 'quarter' was as 
follows: 476 (QO), 1639 (Ql), 4354 (Q2), 4370 (Q3), 4397 
(Q4) and 4633 (Q5); (7) the resulting light curves were fit- 
ted according to their Fourier descriptions and secondary 
trending (usually at the 1-2 sigma level) was removed. 

This method of normalization works well for unmodu- 
lated stars but might explain the observed apparent absence 
of amplitude modulation for V349 Lyr. Owing to its relative 
faintness (V349 Lyr is the faintest star in our sample), and, 
if it is a Blazhko star as suggested by BIO, its long Blazhko 
period and low Blazhko amplitude, we may have to wait an- 
other year or longer to decide if the variations are intrinsic 
or due to instrumental effects. 



2.2.2 Short-cadence data 

In addition to the LC observations, FN Lyr and AW Dra 
were observed at short cadence (SC, 1 minute) in QO and Q5. 
The time baseline for the QO observations was 9.7 d (0.528 d 
to 10.254 d), while the the Q5 observations were made over 
94.7 d (322.98 d to 417.67 d). When combined the total 
number of data points amounted to ~150,000 observations 
per star (14280 from QO, and 136140 from Q5). Processing 
of the data required removal of a small number of outliers 
(cr > 5) and three spurious dips occuring in QO at times 
centered on 3.298 d, 6.261 d and 8.342 d. The SC data were 
detrended in the same way as the LC observations. 



For the eight stars with periods given in the CEOS 
database agreement with the Kepler periods is very good, 
the average difference amounting to only 0.4 x 10~^d. The 
biggest period discrepancies are for V715 Cyg, V784 Cyg 
and AW Dra. The RT differences are discussed in §3.3. Pos- 
sible slowly-changing periods that may not be detectable in 
the Q0-Q5 Kepler data were assessed from O— C diagrams 
and are discussed in §3.4. 

Nine of the Kepler non-Blazhko stars also were observed 
during the course of the ASAS-North survey (see Pigulski 
et al. 2009, Szczygiel et al. 2009, and the ASAS website at 
jhttp:/ /www. astrouw.edu.pl/asas I. For the nine stars ~80 V 
and a similar number of / CCD photometric measurements 
were made over ~500 days in 2006 and 2007. A reanalysis of 
the on-line ASAS data was made and risetimes and total am- 
plitudes, Fourier-fitted (V) magnitudes, and (V — /) colours 
are given in the last four columns of Table 2. If the ASAS- 
South survey is any indication the ASAS-North RR Lyr light 
curves will be "close to the standard V system" (Kovacs 
2005). The ASAS V data are discussed further in §4.5 be- 
low, and the cyclic behaviour of the stars in the HR-diagam 
is discussed in Appendix B. 

Most recently, pulsation periods and preliminary 
Fourier analysis for 16 of the 19 stars were reported by BIO. 
They analyzed the long cadence Kp photometry from QO 
(9.7d interval), Ql (33. 5d) and Q2 (89d) and established 
the non-Blazhko behaviour of most of the stars. The peri- 
ods (typical uncertainties of ±3 x 10^'' d) and Fourier Ai 
amplitudes that they derived are summarized in columns 3 
and 4 of Table 2. 



3 PERIODS AND LIGHT CURVES 

3.1 Previously published periods, risetimes, etc. 

All but one of the 19 stars in Table 1 (the exception being 
KIC 7021124) has been studied previously, either from the 
ground or using the early Kepler photometry. In Table 2 
some of this information is summarized. 

The CEOS RR Lyr database (see LeBorgne et al. 2007) 
at http://rr-lyr.ast.obs-mip.fr/dbrr gives pulsation periods 
for eight of the 19 stars - these are summarized in column 2 
of Table 2. Much of the historical photometry was derived 
from photographic plates taken through various filters and 
on different emulsions, often but not always approximating 
Johnson V and B magnitudes. The database also includes 
HJD times of maximum light, risetimes (RTs), and visual 
and CCD observations acquired over many decades. 

Two historical studies of particular interest are the 
impressive 1953 and 1956 photographic investigations by 
W.J. Miller, S.J. The former study includes an analysis of 
V715 Cyg (then known as Vatican Variable 14, or VV 14), 
while the latter includes analyses of V782 Cyg (VV 22) and 
V784 Cyg ( VV 24) . Both studies were based on photographic 
magnitudes (m-pg) derived from Harvard and Vatican plates 
taken between 1927 (JD 2425148) and 1952 (JD 2434323). 
The number of individual rripg observations were as follows: 
840 for V715 Cyg, 735 for V782 Cyg, and 741 for V784 Cyg. 
Observed times of maximum light are given for V782 Cyg 
(46 epochs) and V784 Cyg (42 epochs) but not for V715 Cyg. 
For all three stars the CEOS catalog adopts the pulsation 
periods and times of maximum light from Miller's investiga- 
tions. 



3.2 Period estimation 

Pulsation periods for the non-Blazhko stars were derived 
mainly from analyses of the Kp data alone. Computa- 
tions were performed with the package PERI0D04 (Lenz & 
Breger 2005), which carries out multifrequency analyses 
with Fourier and least-squares algorithms, and with a ver- 
sion of the CLEAN program (see Nemec, Walker & Jeon 
2009). Details of the derived periods are given in Table 3, 
where the best estimates of the periods and times of maxi- 
mum light are in columns 3 and 4, and column 5 describes 
the specific data that were analyzed. The uncertainties in the 
periods (given in parentheses) are average values obtained 
using the three methods available in PERIQD04, i.e., a theo- 
retically expected uncertainty based on analytical consider- 
ations, a least-squares method, and a Monte Carlo routine 
similar to the Fisher method of randomization described by 
Nemec & Nemec (1985). Typically the periods derived using 
the Q0-Q5 LC data alone are accurate to 1 — 2 x 10~^ d. The 
uncertainty is expected to decrease further as more Kepler 
observations become available. The most uncertain period 
is that for KIC 7021124 which so far has been observed only 
in Ql. Also given in Section (a) of Table 3 are the mean Kp 
magnitudes obtained from the KIC. In column 5 the abbre- 
viations 'LC and 'SC stand for 'long cadence' (30-min) and 
'short cadence' (1-min) Kepler data. Columns 6-11 contain 
Fourier-based results and are described below. 

In general, the new periods compare favourably with 
the CEOS and BIO periods (see Table 2), with the newer 
values being more precise and more accurate. The largest 
differences are for V715 Cyg and AW Dra, both of which 
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Table 2. Other periods, amplitudes, risetimes and ASAS (y)-magnitudes and {V — 7)-colours. 
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have significantly longer periods than given in the GEOS 
catalog. In §3.4 the periods derived directly from the Kp 
data are compared with the values from O— C diagrams. 



3.3 Light curves, amplitudes and risetimes 

Phased light curves for the Kepler non-Blazhko ab-type 
RR Lyr stars are plotted in the left panel of Figure 2, 
assuming the pulsation periods and times of maximum light 
given in Table 3. The light curves have been ordered from 
smallest to largest amplitude (measured from minimum to 
maximum light), with each star offset from the next by 0.5 
mag. This arrangement shows that there is a clear trend in 
light curve shape and risetime (RT), with the lowest ampli- 
tude stars having the longest RTs and the most symmet- 
ric light curves, and the largest amplitude stars having the 
shortest RTs and most asymmetric light curves. 

The largest amplitude star is V368 Lyr, with Atot = 
1.133 mag. This star also has the most asymmetric light 
curve and the shortest period, P = 0.456 d. At the top of 
the stack of light curves the star with the smallest amplitude 
is V782 Cyg with Atot = 0.523 mag; however, its period, 
P = 0.524 d, is not the longest. That distinction belongs 
to three stars with approximately equal periods: AW Dra 
(0.687 d), NR Lyr (0.682 d) and KIC 7030715 (0.684 d). 
These three stars all have intermediate amplitudes, and are 
shown below to have considerably lower metallicities than 
V782 Cyg (see §6.1 and Figs. 7-9). 

The right panel of Figure 2 also shows a comparison 
of the total amplitudes, Atot (Table 1, column 9), and the 
Fourier Ai coefBcients (see next §), both derived from the 
ifp- photometry. Although A\ is approximately proportional 
to Atot it is clear that the relationship is non-linear. When 
a cubic equation was fitted to the graph the resulting stan- 
dard error of the fit amounted to only 3 mmag. It is doubtful 
that this (slight) curvature would have been detected with- 



out the high precision of the Kepler data. In general the 
relationship is roughly approximated by non-linear hydro- 
dynamic models (see §7, in particular the top-left panel of 
Fig. 14). 

The RTs, magnitude ranges, and total amplitudes given 
in Table 1 (columns 7-9) were calculated numerically from 
the fitted Kp light curves. The magnitudes at maximum and 
minimum light are the values at which the slopes of the light 
curves were found to be zero, from which magnitude ranges 
and precise total amplitudes were calculated. The phases at 
the light minima and maxima were used to calculate the 
RTs. 

It is informative to compare the _Kp-based RTs with 
those given in the GEOS database and with RTs based on 
the ASAS ^/-photometry. Since RTs are not given at the 
ASAS on-line website the values given in column 5 of Table 
2 were derived in the same way that the Kp values were com- 
puted. For the nine stars with Kp and ASAS RTs the agree- 
ment is excellent, the mean difference being 0.002 (ASAS 
minus Kp, standard deviation 0.014). For Miller's three Vat- 
ican variables (values are given in parentheses in column 5 
of Table 2) there is reasonable agreement, with the mean 
difference being —0.024 (Miller minus Kp, standard devia- 
tion 0.021). The largest differences occur for NR Lyr and 
V838 Cyg, where the GEOS database gives RT=0.27 and 
0.40, respectively, compared with the averages of the Kp 
and ASAS estimates of 0.147 and 0.132. An independent es- 
timate of 0.141 for NR Lyr comes from our analysis of the 
unpublished V CCD photometry of Benko and Nuspl, which 
agrees well with the Kp value. 



3.4 O C diagrams and Period Change Rates 

O— C diagrams can be used to improve the precision of esti- 
mated pulsation periods, to detect slowly changing periods, 
and to derive accurate period change rates. The longer the 
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Figure 2. (Left) Phased light curves for the 19 non-BIazhko ab-type RR Lyr stars observed with the Kepler telescope. The quantities 
in parentheses following the star names are the total Kp amplitudes (Table 1, column 8), ranging from Atot = 0.52 to 1.13 mag. Almost 
all of the light curves contain over 18000 points with a spread of less than one milli-magnitude. (Right) Graph comparing Atot values 
and the Fourier A\ coefficients (in column 8 of Table 3 below), all derived from the iCp-photometry. The fit to the points is described by 
the cubic equation Ax = 0.443 A^^^ - 0.950 A^^^ + 0.973 Atot - 0.120, with fit standard error a = 0.003 mag. 



time baseline the greater the leverage on the period and the 
greater the probability of detecting curvature signifying a 
changing period. For those Kepler stars with over one year 
of nearly continuous high-precision photometry (Q0-Q5) we 
have seen that it is possible to derive periods accurate to 
~lxlO~^ d using only the Kepler photometry. When the 
Kepler data are combined with the long-baseline historical 
photometry (in some cases as many as 105 years - see for 
example Le Borgne et al. 2007) even more accurate periods 
(and possible dP/dt values) can be calculated. 

Figure 3 shows O— C diagrams for the eight Kepler 
non-Blazhko stars that have ephemeris information in the 
GEOS database. The ordinate represents the "observed mi- 
nus calculated" phase at which maximum light occurs, as- 
suming the GEOS period and time of maximum light given 
in parentheses at the top of each panel; the abscissa is the 
Heliocentric Julian Date0 At each epoch the O— C value 
was calculated by taking the observed time of maximum 
light minus the assumed GEOS time of maximum light (to), 
and dividing the resultant by the assumed (constant) GEOS 
period (P). This gave the number of pulsation cycles since 
to, and the remainder after subtracting an integer number 
of cycles (E) is the observed fractional phase shift from the 
expected value. For five of the stars linear fits to the obser- 
vations are given, a positive (negative) slope indicating that 
the assumed period is too short (long). Quadratic fits are 
shown for V782 Cyg, FN Lyr and AW Dra (shown twice, 
the second graph representing the last ~13 years), where 

^ Since the offset between Barycentric Julian Dates (BJDs) used 
by Kepler are very small compared to the HJDs usually used by 
ground-based observers, the two types of dates have been used 
interchangeably in this paper. 



the upwards curvature indicates an increasing period, the 
period change rate, dP/dt, being given by 2Pc^ where P is 
the assumed GEOS period, and c is the curvature coefficient 
in the quadratic equation. None of the stars represented in 
Fig. 3 was studied by Le Borgne et al. (2007). 

For the three Vatican Variables (V715 Cyg, V782 Cyg, 
V784 Cyg) there is an ~60 year gap between the 2009-10 
Kepler data and the earlier observations. Over such a long 
time interval one cannot hope to keep track of cycle counts 
and there is little hope of using the O— C diagrams to ei- 
ther improve the period or calculate dP/dt. To illustrate the 
cycle-count problem the panels for V715 Cyg and V784 Cyg 
show two (of many) possible linear solutions, where, in all 
the cases there is a large O— C shift between the Kepler and 
earlier epochs0 However, comparison of the assumed periods 
with the periods derived from the Kepler data do give some 
indication of either a period change or incorrect period. For 
V715 Cyg, the ifep/er period, 0.4707045(1) d, is significantly 
longer than the GEOS period, 0.47067298(43) d, suggest- 
ing a period increase. And for V784 Cyg the Kepler period, 
0.5340941(1) d, is significantly shorter than that derived by 
Miller (1956), 0.5341026(11) d. 

Miller's (1956) data for V782 Cyg provides 46 times 
of maximum light at epochs between 1928 and 1954, and 
with these Miller derived P = 0.52363383 d (adopted in the 
GEOS database). Our analysis of the ~18419 Q1-Q5 Kp 



3 Miller (1953) did not bin his V715 Cyg (VV 14) observations 
into separate epochs, and the GEOS database simply adopted for 
its zero-phase epoch the time of maximum light given in the Miller 
ephemeris. Future period change rate studies of this star might 
want to bin Miller's photometry into separate epochs rather than 
just let one epoch (1948) represent the extended 1927-1952 data. 
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Figure 3. O— C diagrams for the eight Kepler non-Blazhko stars with previous photometry in the GEOS atlas (AW Dra is shown twice, 
with different time seales). In each panel the O— C value derived from the Kepler data is plotted as a solid black square at HJD ~55000, 
and the crosses are the O— C values at the various previous epochs. Error bars for the Kepler points are comparable to the size of the 
square and have not been plotted. The fitted lines (or curves) plotted in each panel, the equations of which are given in each graph, are 
discussed in the text. For V715 Cyg and V784 Cyg two, of many, possible solutions are plotted. 



photometry gave P = 0.5236377(1) d, which is 3.9xl0~'^d 
longer than the Miller period (see column 3 of Table 2). 
Using the time of maximum light given in column 5 of Ta- 
ble 1 the O— C phase shift of ~0.20 (nearest) is undoubtedly 
significant. The quadratic solution (shown in Fig. 3) gives 
dP/dt = 0.06 d/Myr, which predicts, for 2010, the period 
0.5236379 d, a value that is only 0.2xl0"'^d longer than the 
measured /("epZer period. While this agreement strongly sup- 
ports the increasing period it is also possible that the period 
changed abruptly in the 60-year gap. 

For V838 Cyg a total of 31 times of maximum light are 
available in the GEOS database: 26 epochs betwen 1958 and 
1968, a single epoch in 1977, and four recent epochs from 
2008-09. Although the 10- and 30-year gaps may be prob- 
lematic, the Kepler ephemeris gives an O-C value near zero 
for every epoch, suggesting either a constant (or very slowly 
changing) period since 1958. A linear fit to the O— C dia- 
gram (Fig. 3, bottom left panel) predicts a (constant) true 
period P = 0.4802797(2) d, which matches, to within the 
uncertainties, the Kepler-hased period P = 0.4802799(1) d 
and the GEOS period P = 0.4802795 d. A parabohc fit sug- 



gests dP/dt — +0.05±0.04 d/Myr, with the assumed period 
0.4802799 d occurring at BJD 48438. The fit standard error 
in both cases is only 0.021 (phase units). 

For NR Lyr (Fig. 3, top left panel) the 31 historical 
points and the new Kepler data are consistent with a shorter 
period than assumed in the GEOS database. The negative 
slope suggests a period shorter by 0.46xl0~^d than the as- 
sumed P(GEOS), i.e., P(true) = 0.6820288 d; this period 
is significantly longer than the estimate based on the 18358 
Kepler Q1-Q5 data points, 0.6820264 (2) d. 

The 13 historical points and the Kepler data for NQ Lyr 
produce consistent period estimates, with a difference of 
only O.lxlO^^d between the two periods. Thus the as- 
sumed (GEOS) period of 0.5877888 d is to be compared 
with 0.5877887(1) d derived using the 18759 Kepler Q0-Q5 
data points. 

The O— C diagram for FN Lyr (middle right panel of 
Fig. 3) assumes P = 0.52739716 d and to = 33829.402 
(GEOS values), and includes 34 O— C values from as early 
as 1907 (t = 17793) to as recent as 2009 {t = 54995). AU 
but one of the phase shifts is near zero (a visual observa- 
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tion from 1937 was considered an outlier and excluded from 
our analysis). The diagram includes a point from the 1966 
photoelectric ^/-photometry of Bookmeyer et al. (1977) at 
t = 39416.6457, and a point from unpublished 2006-08 CCD 
F-photometry by Layden at t = 54288.7454 (see §4.4 for fur- 
ther discussion of these data) . These two points, which differ 
in time by ~40 years, have been plotted as open circles and 
differ in phase by ~0.055. The O— C value derived for FN Lyr 
from the Kepler data (epoch 2009) is almost identical to the 
Layden O— C value (see Table 2). Considering that the to- 
tal baseline is longer than 100 years the O— C phase shifts 
are very near zero. Nevertheless there appears to be a slight 
upward trend, and possibly a small amount of curvature. A 
parabolic solution gives dP/dt=0.004±0.004 d/Myr, with 
the GEOS period occuring at t=18100, a residual standard 
error of the fit of 0.011, and a predicted period 0.5273975 d 
for the current epoch (2010-11). Since the curvature term 
does not differ significantly from zero (p-value= 0.27) the 
question of a changing period remains unanswered. A lin- 
ear solution is also consistent with the data, in which case 
a revised constant period of 0.52739741(5) d is suggested 
by the data. Again, the residual standard error of the fit 
was only 0.011. When the 136140 SC observations made in 
Q5 were combined with the 14280 SC observations made in 
QO, a period 0.52739845(1) d with residual standard error 
0.60 mmag, was derived, borderline consistent with both the 
constant period and small AP /At hypotheses. 

Finally, the GEOS data for AW Dra have many epochs 
with O— C phases that appear to be outliers. These include 
a single point at t ~20000 with O— C~ —0.20, and a clump 
of ~20 points at t ~52000 with 0-C~ -0.30. If these val- 
ues are correct then they present a challenge for modelling 
period variations. We have chosen to discard them as out- 
liers leaving the O— C diagram shown in the bottom mid- 
dle panel of Fig. 3. The parabolic fit suggests an increas- 
ing period amounting to 0.32 d/Myr, and predicts a period 
0.687214 d for the Kepler Q5 epoch. However, close inspec- 
tion of the diagram shows that the last three points, one of 
which is that derived from the Kepler data, deviate substan- 
tially from the fitted curve. If the analysis is restricted to 
epochs since HJD 49000 (the first point is from the Castel- 
lani et al. 1998 study of AW Dra) and if the suspected out- 
liers identified above are excluded, then the O— C diagram is 
as shown in the bottom right panel of Fig. 3. In this case, the 
period increases smoothly at the rate AP /At — 3.76 d/Myr, 
with the assumed period of 0.6871941 d having occurred at 
HJD 49267. For the Kepler Ql and Q5 LC data (6117 points) 
the predicted period is 0.6872156(2) d. And the period from 
the combined QO and Q5 SC data (149222 points) is similar, 
0.6872162(1) d, with a fit standard error 7.3 mmag. Both pe- 
riods agree, are significantly longer than the GEOS period, 
and are consistent with AW Dra having an increasing pe- 
riod. Comparing this period with the predicted periods the 
new data favour the smaller AP/At value. 

To summarize this section, periods based on O— C di- 
agrams were calculated for NR Lyr and NQ Lyr, both of 
which seem to have constant periods, and AP/At values 
have been calculated for FN Lyr and AW Dra. AW Dra 
is the only star for which we can at present be confident 
that the period is changing (assuming that the outliers re- 
ally should be omitted from consideration) - depending on 
the assumed data set the rate is either dP/dt=3.786 d/Myr 



(for the more recent data), or 0.32 d/Myr (all data) with 
the possibility of variations on top of this. For FN Lyr the 
period either is very nearly constant or is increasing at a 
very slow rate. For the other four stars with GEOS times of 
maximum light the large gaps in the observation record pre- 
vent improvement over the GEOS or Kepler periods. The 
ongoing high-precision Kepler observations should remedy 
this lack of knowledge (for the modern era), not only for the 
eight Kepler non-Blazhko stars in the GEOS database but 
also for all the Kepler non-Blazhko stars. 



4 FOURIER ANALYSIS 

The use of Fourier methods to characterize the shapes of 
RR Lyr light curves began with the work of Simon and his 
collaborators (Simon & Lee 1981, Simon & Teays 1982, Si- 
mon & Clement 1993) and has continued through the more 
recent papers by Kovacs, Jurcsik, Morgan and others. The 
goal of these studies has been to establish empirical correla- 
tions employing a minimum number of observed parameters 
(periods, amplitude ratios, phase parameters) from which 
physical characteristics can be derived (see §6). Calibra- 
tion of the equations follows from independent observations 
(high dispersion spectroscopy, parallax measurements, etc.) 
and from theory {e.g., hydrodynamic pulsation models, stel- 
lar evolution models). After Simon (1988) discovered that 
the Fourier phase parameter <^2i decreases with decreasing 
metallicity (for field RRab stars with P < 0.575d), Kovacs & 
Zsoldos (1995) proposed that correlations involving period 
and the Fourier phase parameter (^31 can be used to derive 
[Fe/H] values. The first application of this method to field 
and globular cluster RR Lyr stars was the study by Jurcsik & 
Kovacs (1995,1996), and subsequent calibrations have been 
made by Jurcsik (1998, hereafter J98) and Kovacs & Walker 
(1999, 2001). From the beginning Jurcsik & Kovacs (1995) 
warned that the Fourier (^31 method "is not applicable for 
the estimation of [Fe/H] in peculiar stars (e.g., in Blazhko 
variables, highly evolved stars ...)", and subsequently JK96 
and Kovacs & Kanbur (1998) introduced a 'compatibility 
test' for identifying 'peculiar' stars. Indeed, when such stars 
were omitted from the ASAS-South database, agreement 
between Fourier-based metallicities and low-dispersion spec- 
troscopic metallicities reached the 0.16 dex level (see Kovacs 
2005). 

4.1 Coefficients and fitted light curves 

Fourier decomposition of the Kepler light curves was per- 
formed for each of the 19 non-Blazhko stars by fitting the 
following Fourier sine series to the observed photometry: 

F 

m{t) ^ Ao + ^AiSm[iujQ{t ~ to) + (t>i\, (1) 

i = l 

where m(t) is the apparent magnitude (either Kp for the 
Kepler Aa±a., or mv(t) for the ground based V-photometry), 
F is the number of fitted terms, luq is the angular pulsation 
frequency of the star (=27r/o, where fo~Po^^), t is the ob- 
served time of the observation (BJD-54953 for the Kepler 
data, HJD for the ground-based V photometry), to is the 
time of maximum light (used to phase the light curves so 
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that maximum light occurs at zero phase), and the Ai and 
(^i are the amplitude and phase coefficients for the individ- 
ual Fourier terms. The Fourier calculations were made with 
two FORTRAN programs, one kindly provided by Dr. Geza 
Kovacs, and the other by Dr. Pawel Moskalik. The assumed 
pulsation periods were those derived either from the Kepler 
data or from the period change rate analysis (§3.4), and the 
to values were calculated numerically. The derived periods 
and to values are summarized in Table 1 (columns 4,5) and 
the details given in Table 3. 



4.2 Kepler ifp-photometry 

Initially the Kepler magnitudes were fitted with a Fourier se- 
ries having 7-15 terms, the larger number being needed when 
the skewness is greatest, i.e., the risetime is shortest. Inspec- 
tion of the residuals from the fitted light curves revealed 
that because most of the stars show a sharp rise to max- 
imum light, with additional detailed bumps and features, 
many more terms were needed. Figure 4 illustrates 38-term 
Fourier fits for two non-Blazhko stars, AW Dra on the left 
and FN Lyr on the right. These two stars have both Kp- 
photometry (top diagrams) and high-precision ground-based 
V-band photometry (bottom diagrams). For AW Dra the 
Kepler data are Q5 SC observations, i.e., brightness mea- 
surements every minute over 94.7 days, resulting in 135380 
data points; and for FN Lyr the Kepler data are Q1-Q5 LC 
observations, i.e., a measurement every 30 minutes over 417 
days, resulting in 18338 points. In both cases the standard 
deviation of the fit to the Kepler data is less than one milli- 
magnitude (see column 6 of Table 3). However, the residuals 
still show systematic variations with the largest residuals oc- 
curing on the rise to maximum light where the slope of the 
light curve is steepest. To achieve residuals showing nor- 
mally distributed white noise perhaps 50-100 terms would 
be needed. 

For each of the program stars epoch-independent phase 
differences c^fi = (f>l — i(f>\, and amplitude ratios Rn — 
Ak /Al , were computed from the Fourier coefficients derived 
from the _Kp-photometry. The superscripts "s" and "c" sig- 
nify phases and phase-parameters computed with sine and 
cosine series, respectively {e.g. Fig. 7 below). The sine and 
cosine (f>3i phase parameters differ by tt radians, that is, 
031 — 031 ~ 3.14159, and the (f)2i parameters differ by tt/2 
radians, i.e., 021 = ^ii + 1.5708. 

Table 3 also contains the Fourier parameters needed for 
computing physical characteristics. In Section (a) results de- 
rived from various combinations of short- and long-cadence 
Kp data (see column 5) are given. The Ai coefficients (col- 
umn 7) are practically identical to the Ai coefficients cal- 
culated by BIO using the Q0-Q2 photometry (given in col- 
umn 6 of Table 2) . Columns 8 and 9 contain the amplitude 
ratios, R21 and 7i3i; and columns 10 and 11 contain the 
Fourier phase parameters, (jiii ^-nd (ji^i. The accuracy of the 
ii'p-based Fourier parameters is at least one part in 1000, 
and higher in many cases. For example, the parameters de- 
rived for FN Lyr from a 38-term fit of the 18338 Q1-Q5 LC 
data points have formal values (including uncertainties) as 
follows: R21 = 0.44208(2), Tisi =0.34610(2), 0^1=2.32417(7), 
and 031=4.81916(9). Correlations and graphical representa- 
tions of these parameters are discussed below. 



4.3 CycIe-to-CycIe variations and stationarity 

To examine cycle-to-cycle variations in the Kepler light 
curves the Fourier parameters were re-calculated after bin- 
ning the data so that each bin included data for one pul- 
sation period {i.e., each bin represented a single pulsation 
cycle). For the long cadence (30 min) data there were ~24 
observations per bin for an RRab star with a period of 12 
hours, and fewer (more) observations per cycle for stars with 
shorter (longer) periods. None of the 19 stars exhibits the re- 
cently discovered "period doubling" effect, thus confirming 
the result found earlier by Szabo et al. (2010). 

For the Fourier calculations 'direct Fourier fitting' (dff ) 
rather than 'template Fourier fitting' (tff) methods were 
used (see Kovacs & Kupi 2007), and for each star the re- 
sulting time series were plotted for four Fourier parameters: 
Al, the first term in the Fourier series (see right panel of 
Fig. 2); 01, the phase of the first term; R21, the ampli- 
tude ratio A2/A1; and 03i, the Fourier parameter found 
by Simon, Kovacs and others to be one of the most sig- 
nificant variables for deriving physical characteristics. Typ- 
ical time series are illustrated in Figure 5, on the left 
for the unmodulated RRab star NQ Lyr having intermedi- 
ate amplitude (Ai=0.279 mag) and intermediate brightness 
{{Kp} = 13.075 mag), and on the right for V783 Cyg, a 
low-amplitude Blazhko star with the shortest Blazhko pe- 
riod {Pb = 27.7 ± 0.04 d - see BIO) among the stars in our 
sample. 

For NQ Lyr linear trend lines were fit to each time series 
and in almost every case the slope is zero to within the 
systematic and random uncertainties. The mean values for 
the four parameters are as follows: 03i = 5.0958 ± 0.0023, 
with residual standard error, a = 0.002; R21 = 0.4710 ± 
0.0006, with a = 0.0006; 0i = 3.961±0.001, with a = 0.001; 
and Al = 0.28016 ± 0.00014 mag, with a = 0.0001. Since 
the other non-Blazhko stars show random variations of the 
Fourier parameters similar to those shown here for NQ Lyr 
these means and errors provide a measure of the typical 
uncertainties. 

Inspection of the NQ Lyr time series in Fig. 5 shows that 
the stability over the ~420 d interval (Q0-Q5) is exceptional, 
not only for each Fourier parameter but for the ensemble of 
parameters. To better illustrate the remarkable stationarity 
of the light curves for all 19 sample stars a set of 'animated 
gif light curves has been prepared and these are available 
in the electronic version of the Journal. 

For the Blazhko star V783 Cyg the time plots of the four 
Fourier parameters are in striking contrast to the stochas- 
tic noise exhibited by NQ Lyr and by the candidate non- 
Blazhko star V349 Lyr. For V783 Cyg pronounced periodic 
variations are seen over the 27.7 d Blazhko cycles. There 
does appear to be a slight discontinuity at around 140 d, 
but this is certainly a data processing defect. Over the 11 
Blazhko cycles seen there is very little difference from cycle- 
to-cycle; this is not always the case for Blazhko stars where 
large cycle-to-cycle differences often are seen. 



4.4 High-precision V-photometry 

By comparing the light curves and Fourier parameters for 
the Kp and V photometric systems it is possible to test (in a 
limited way) the hypothesis that the two broad-band filters 
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Figure 4. Light curves and residual plots for AW Dra (left) and FN Lyr (right). The top diagrams show Kp photometry and the 
bottom diagrams V^-band photometry, and within each diagram the light curves are on top and the residual plots on the bottom. The 
Kp photometry is SC(Q5) data for AW Dra and LC(Q1-Q5) data for FN Lyr. The x-axis labels give the assumed pulsation periods (Po) 
and times of maximum light (to)- The diagram labels indicate the data plotted and the number of terms in the fitted Fourier series, and 
the number of data points (N) and standard deviations of the fits (cr) arc stated in each diagram. 



give similar results. If this hypothesis is true then the Fourier 
parameter correlations established for the V-band system by 
Simon, Kovacs, Jurcsik, Sandage, and others might be ap- 
plied to the Kp mags and be used to derive first-order physi- 
cal characteristics (see §6). Unfortunately there is very little 
high-precision published ground-based V-band photometry 
for our Kepler stars. Three exceptions are NR Lyr, AW Dra 
and FN Lyr. Using the available V-photometry for these 
stars Fourier parameters were derived and compared with 
the Kp parameters. The results are given in the middle sec- 
tion of Table 3 (Section b). 



4.4.1 AW Dra 

The Ught curve and residuals derived from the V- 
photometry of Castellani et al. (1998) are plotted in the 
lower left diagram of Figure 4. The 112 points are from CCD 
observations made in 1995 with the 72-cm Teramo/Scuola 
Normale Telescope, modelled with a 7-term Fourier series fit. 
The shape of the V light curve is very similar to the light 
curve based on the Kepler data, and the Fourier parame- 
ters (Table 3) support this observation. Notice also that, 
although the V-data have a residual standard error of only 
0.0116 mag, the residuals are ~ 16x larger than the Kp- 
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Table 3. Pulsation and Fourier parameters for tlic Kepler non-Blazhko ab-type RR Lyr stars 



Star (Kp) Period to Observational cr Ai R21 -R31 'j>2i '/'ii 

[mag] [day] [BJD] data (No.pts) [mmag] [mag] [rad] [rad] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



(a) Results from analysis of the Q0-Q5 Kp photometry 



NR Lyr 


12.683 


6820264('2"1 


54964 7381 


LC-01-05 fl8333"l 


0.69 


0.266 


0.456 


0.352 


2.416 


5.115 


V71 'i Cvp- 


16.265 


n 47n7n494('4"i 


54Qfi4 




1.74 


0.338 


0.479 


0.358 


2.314 


4.901 


V782 CvET 


15.392 


5236377('l ") 


54Qfi4 5059 


LC 01-05 flSSSTl 


0.79 


0.190 


0.488 


0.279 


2.777 


5.808 


V784 Cyg 


15.370 


0.5340941(1) 


54964.8067 


LC:Q1-Q5 (18364) 


0.96 


0.234 


0.487 


0.253 


2.904 


6.084 


KIC 6100702 


13.458 


0.4881457(2) 


54953.8399 


LC:Q0-Q4 (14404) 


0.66 


0.209 


0.493 


0.279 


2.743 


5.747 


NQ Lyr 


13.075 


0.5877887(1) 


54954.0702 


LC:Q0-Q5 (18759) 


0.65 


0.280 


0.471 


0.356 


2.389 


5.096 


FN Lyr 


12.876 


0.52739845(1) 


54953.2690 


SC:Q0-I-Q5 (149925) 


0.60 


0.380 


0.445 


0.354 


2.321 


4.817 




12.876 


0.527398471(4) 


1 54953.2690 


LC:Q1-Q5 (18338) 


0.68 


0.379 


0.442 


0.346 


2.324 


4.819 


KIC 7021124 


13.550 


0.6224926(7) 


54965.6471 


LC:Q1 (1595) 


1.10 


0.283 


0.512 


0.351 


2.372 


5.060 


KIC 7030715 


13.452 


0.6836137(2) 


54953.8434 


LC:Q0-Q5 (18802) 


0.71 


0.231 


0.494 


0.303 


2.683 


5.606 


V349 Lyr 


17.433 


0.5070740(2) 


54964.9555 


LC:Q1-Q5 (18314) 


3.24 


0.346 


0.450 


0.352 


2.328 


4.845 


V368 Lyr 


16.002 


0.4564851(1) 


54964.7828 


LC:Q1-Q5 (18273) 


1.57 


0.405 


0.464 


0.341 


2.272 


4.784 


V1510 Cyg 


14.494 


0.5811436(1) 


54964.6695 


LC:Q1-Q5 (18394) 


0.82 


0.345 


0.473 


0.355 


2.389 


5.068 


V346 Lyr 


16.421 


0.5768281(1) 


54964.9211 


LC:Q1-Q5 (18362) 


2.83 


0.330 


0.473 


0.352 


2.372 


5.060 


V350 Lyr 


15.696 


0.5942369(1) 


54964.7795 


LC:Q1-Q5 (18326) 


1.67 


0.340 


0.485 


0.342 


2.389 


5.124 


V894 Cyg 


13.293 


0.5713866(2) 


54953.5627 


LC:Q1-Q5 (18362) 


0.91 


0.377 


0.490 


0.338 


2.364 


5.067 


V2470 Cyg 


13.300 


0.5485894(1) 


54953.7808 


LC:Q0-Q5 (18794) 


0.79 


0.220 


0.488 


0.282 


2.745 


5.737 


V1107 Cyg 


15.648 


0.5657781(1) 


54964.7532 


LC:Q1-Q5 (18373) 


0.99 


0.280 


0.495 


0.350 


2.421 


5.196 


V838 Cyg 


13.770 


0.4802799(1) 


54964.5731 


LC:Q1-Q5 (18241) 


1.22 


0.393 


0.465 


0.349 


2.300 


4.853 


AW Dra 


13.057 


0.6872158(6) 


54954.2160 


SC:QO (14240) 


0.55 


0.305 


0.527 


0.348 


2.731 


5.563 




13.053 


0.687217(1) 


54954.2160 


LC:Q1 (1614) 


0.69 


0.306 


0.524 


0.342 


2.730 


5.561 




13.053 


0.6872158(2) 


54954.2160 


LC:Q5 (4474) 


0.79 


0.306 


0.524 


0.343 


2.728 


5.557 




13.053 


0.68721632(3) 


54954.2160 


SC:Q5 (135380) 


0.71 


0.308 


0.527 


0.347 


2.729 


5.558 






(b) Results 


from (re-)analysis of high-precision ground-based V photometry 








NR Lyr 




0.6820264 


54964.7381 


Benko & Nuspl (183) 


5.7 


0.322 


0.445 


0.343 


2.32 


4.94 


FN Lyr 


12.810 


0.52739716 


54953.2630 


Layden unpub. (124) 


30.1 


0.436 


0.433 


0.344 


2.25 


4.67 


AW Dra 


12.823 


0.6871941 


49918.2925 


Castellani et al. (112) 


11.6 


0.357 


0.510 


0.353 


2.63 


5.43 




(c) Results from rc- 


-analysis of the ASAS V photometry l|http://www 


.astrouw.edu. pl/asas| 






NR Lyr 


12.44 


0.68200(5) 


54964.7383 


190827-1-3848.8 (82) 


61 


0.24 


0.45 


0.35 


2.2 


4.9 


KIC 6100702 


13.64 


0.48817(2) 


54953.8398 


185038-1-4125.4 (79) 


71 


0.23 


0.59 


0.32 


2.4 


4.8 


NQ Lyr 


13.36 


0.58779(3) 


54954.0702 


190749-1-4217.9 (68) 


68 


0.30 


0.56 


0.40 


2.4 


5.2 


FN Lyr 


12.78 


0.52740(3) 


54953.2690 


191022+4227.6 (82) 


76 


0.41 


0.49 


0.36 


2.3 


4.8 


KIC 7030715 


13.24 


0.68362(3) 


54953.8434 


192325-1-4231.7 (84) 


103 


0.26 


0.56 


0.35 


2.4 


5.4 


V894 Cyg 


12.92 


0.57140(2) 


54953.5625 


193301-1-4614.3 (76) 


108 


0.37 


0.41 


0.38 


2.2 


4.8 


V2470 Cyg 


13.53 


0.54857(1) 


54953.7808 


191958-1-4653.3 (78) 


62 


0.25 


0.42 


0.33 


2.7 


5.6 


V838 Cyg 


14.25 


0.48030(2) 


54964.5731 


191404-1-4812.1 (68) 


105 


0.44 


0.39 


0.38 


2.2 


4.7 


AW Dra 


12.85 


0.68721(1) 


54954.2158 


190048-1-5005.5 (86) 


49 


0.37 


0.53 


0.36 


2.5 


5.4 



residuals (Figure 4, upper left panel). On the other hand, 
the risetimes are quite similar: 0.165 for the /fp-photometry 
and 0.174 for the ^/-photometry. The joined points in the RT 
graphs below (lower left panel of Figure 9, and left panel of 
Figure 10) give a visual representation of these differences. 

i.4.2 FN Lyr 

For FN Lyr, the V photometry shown in the lower right pan- 
els of Figure 4 are from two sources: Bookmeyer et al. (1977) 
and Layden (2010, unpublished). The Bookmeyer et al. data 
comprise 71 UBV photoelectric observations made in 1966 
(21) and 1969 (50), and the Layden T/-photometry consti- 
tute 124 VR CCD observations made in 2006-08. When the 
two data sets, which are separated by ~40 years, were com- 
bined and plotted as single hght curve [i.e., the same pe- 



riod and time of maximum light) a phase shift amounting 
to 5.5% of the period (0.029 d) was seen, with the Lay- 
den data having the larger phases. The shift corresponds 
to the offset of the two O— C values seen in Fig. 3 (mid- 
dle right panel, where the O— C values for the two data 
sets are plotted as open circles and are labeled '1966-69' 
and '2006-08'). Our period change rate analysis suggested 
that the shift is probably due either to the assumed period 
(0.52739716 d) being too small or to a slowly increasing pe- 
riod {dP/At= 0.004 ± 0.004 d/Myr). 

Figure 4 (lower right) shows the composite V light curve 
for FN Lyr, where the Layden data (solid black squares) 
are plotted alongside 'the Bookmeyer data shifted in phase 
by +0.055' (red open circles). There is a reasonably good 
match between the Layden and Bookmeyer light curves, 
the chief differences being: (1) only two (very close) Lay- 
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NQ Lyr (KIC 6V63132, P=0.5B77BB7d. t0=64953+ 1 ,070B) 
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Figure 5. Variations in time of tiie Fourier parameters ^gj, R21, 0f and Ai (all sine-series) derived for tlie typical non-Blaziiko star 
NQ Lyr (left panel) and for the known low-amplitude Blazhko star V783 Cyg (right panel). Eacii point corresponds to the variable of 
interest as determined from a single pulsation cycle. 



den data points were acquired near maximum light, and they 
are fainter than the Bookmeyer observations near maximum 
hght; (2) the slope of the rise to maximum light is shallower 
for the Layden data than for the Bookmeyer data; and (3) 
the Layden data are uniformly distributed in phase while 
the Bookmeyer data show two large gaps on the descent to 
minimum light. 

Fourier coefficients and parameters were calculated us- 
ing the composite 'Layden plus shifted-Bookmeyer' data, 
and using the 124 Layden V data points only (the for- 
mer were found to be less reliable than the latter). Owing 
to the large gaps it was not possible to obtain a reliable 
Fourier fit using only the Bookmeyer data. Table 3 gives 
the Fourier decomposition results based on the Layden-only 
y-photometry. The l/-residuals shown in the bottom-right 
panel of Figure 4 were calculated with respect to the Layden- 
only fit. The Layden data have a residual standard error of 
0.03 mag, while the fit for the Bookmeyer data is poorer 
(~0.06 mag). In both cases (as for the Kepler data) the 
largest residuals occur on the rise to maximum light. The 
observed smaller total amplitude and longer risetime for the 
Layden V photometry compared with the Bookmeyer pho- 
tometry might be expected if FN Lyr is a low-amplitude 
Blazhko variable. The available data are insufficient to check 
this possibility but as additional Kepler observations are ac- 
quired the answer may be evident. 



4.4.3 NR Lyr 

In 2008 B, V, R, I CCD observations of NR Lyr were made 
by J. Benko and J. Nuspl using the Konkoly telescope. At 
present the 183 data points per filter are on the instrumen- 
tal system of the telescope; however, the photometric colour 
constant is very small and the V data differ from the John- 
son B, V system by only a zero point shift. Because the 



data completely cover the phase range and are of high pre- 
cision the V data were Fourier analyzed (using the period 
and to-value favoured by the Kepler data) and the results 
summarized in Table 3. When the Fourier parameters are 
compared with those for FN Lyr and AW Dra (see Figures 
6, 9 and 10) all three sets of observed offsets are in excellent 
agreement (see §5.2 below). 



4.5 ASAS V'-photometry 

For the nine non-Blazhko RR Lyr stars in the ASAS-North 
survey (see §3.1) the on-line ^-photometry was re-analyzed 
and the results reported in section (c) of Table 3. The de- 
rived {V} magnitudes (column 2) have uncertainties of ±0.01 
and range from 12.44 to 14.24. A comparison of the (V^) 
and (Kp) mags shows the following trend over the observed 
range: (V") = (1.45±0.24)(A'p) - (5.97±3.20). Untfi further 
V photometry fainter than 14"^ mag is obtained it is unclear 
whether this trend of V mags fainter than the Kp mags will 
continue. Pulsation periods (with uncertainties in parenthe- 
ses) derived from our re-analysis of the V data are given in 
column (3); in general these agree with the periods given at 
the ASAS website and with those derived from the Kepler 
data (see Tables 2-3). Columns 4-6 contain, respectively, the 
assumed time of maximum light, the ASAS name of the star, 
and the spread of the Fourier fit about the mean light curve 
(found to range from 49 to 108 mmag). Fourier parameters 
derived using the on-line V data are given on the right side 
of Table 3. These were computed using the accurate Kepler 
periods and times of maximum light. The uncertainties typ- 
ically are ±0.04 for R21, ±0.01 for .R31, ±0.08 for (/)2i, and 
±0.1 for 031, and in several cases are larger. 
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0.54 




Figure 6. Four panels comparing various Fourier parameters derived from the ifp-photometry for the 19 Kepler non-Blazhko stars. 
The open diamonds for AW Dra, FN Lyr and NR Lyr were derived from the V-photometry of Castellani et al, (1998), Layden (2010, 
unpublished) and Benko & Nuspl (2010, unpublished), respectively. 



5 FOURIER PARAMETER CORRELATIONS 



5.1 Kp Correlations 

Figure 6 shows four graphs relating the Fourier parame- 
ters Al, R31, R21, 021 and 4>3i, derived from the sine-series 
decomposition of the ifp-photometry of the 19 Kepler non- 
Blazhko stars. Each star is represented by a solid black 
square, and is labelled for easy identification. Also plotted 
in Fig. 6 are the F-band Fourier parameters for AW Dra, 
FN Lyr and NR Lyr (§4.3) which are shown as open dia- 
monds connected to the corresponding Kp points by lines. 
The offsets between the Kp and V points appear to be sys- 
tematic (i.e., generally go in one direction) and are not very 
large, even in the case of Rsi which is the least certain of the 
parameters. Perhaps this relationship between V and Kp is 
not surprising because both passbands are quite wide. 

The upper-left panel of Fig. 6 shows that there is very 
little correlation between R21 and period (except possibly 
a slight increase in R21 with increasing period). The R21 
values from the V photometry are all slightly smaller than 



from the Kp photometry, and the R21 values for FN Lyr and 
NR Lyr are significantly smaller than that for AW Dra. 

The upper-right panel shows 14 stars (including 
AW Dra, FN Lyr and NR Lyr) with R31 between 0.34 and 
0.36, and five stars with R^i values smaller than 0.31. In par- 
ticular, AW Dra, FN Lyr and NR Lyr are seen to have sim- 
ilar i?3i values even though the R21 values for FN Lyr and 
NR Lyr are significantly smaller than that for AW Dra. Al- 
though the overall bimodal distribution is striking, we shall 
see below (upper right panels of Figs. 7 and 8) that RR Lyr 
stars with relatively low R31 values are quite common and 
can be simulated with hydrodynamic models (§7). 

In the lower-left panel the stars separate into high- and 
low-(/)2i stars, with a discontinuity at j4i~0.25 mag. AW Dra 
and V784 Cyg (suspected of having a relatively high metal- 
licity) appear to stand out from the other stars; even when 
compared with RR Lyr stars in globular clusters (§5.3) these 
two stars appear extreme. 

Finally, the lower-right panel shows that 02 1 is roughly 
proportional to 03i, where the linear relationship is given by 
(1)21 = 0.509 031 - 0.181. For AW Dra, FN Lyr and NR Lyr 
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Figure 7. Four panels comparing the Fourier parameters (transformed to tiie V-band) for tlie Kepler non-BIazliko stars (large black 
squares) and for the 177 RR Lyr stars located in several galactic and LMC globular clusters (small black dots). The cluster RR Lyr stars 
(from Kovacs & Walker 2001) are in the globular clusters M2 (12 stars), M4 (4), M5 (12), M55 (4), M68 (5), M92 (5), NGC 1851 (11), 
NGC 5466 (6), NGC 6362 (12), NGC 6981 (20), IC 4499 (49), Ruprecht 106 (12), NGC 1466 (8), Reticulum (8) and NGC 1841 (9). 
The points for the 19 most metal-poor globular cluster stars are surrounded by blue squares, and the points for the 39 most metal-rich 
globular cluster stars are circled in red. 



the phase parameters are smaller for the V data than the 
Kp data, and the offsets are parallel to the fitted line. 

5.2 Kp-V offsets 

Since the offsets between the Kp and V Fourier parame- 
ters for AW Dra, FN Lyr and NR Lyr appear to be small 
and systematic (Figs. 6,9,10; Table 3), transformation from 
Kp to V should be possible using the high-precision V pa- 
rameters (inclusion of the ASAS parameters would not have 
been helpful owing to their relatively large uncertainties). 
Transformation to the V system allows derivation of physical 
characteristics through the application of well-established 
V-band Fourier relations (see §6 below). Because the three 
stars exhibit a range of periods (0.53 d to 0.69 d) and hght 



curve shapes (^li = 4.8 to 5.6) and mean colours (see Fig. Bl 
below), confounding effects that depend, for example, on lo- 
cation in the 'instability strip', are likely to have been re- 
vealed. On the other hand, if the offsets depend on [Fe/H] 
such an effect might not have been detected because all three 
stars appear to have low metal abundances (see column 3 of 
Table 4). 



Our approach to estimating the Kp-V offsets was simply 
to average the observed differences for AW Dra, FN Lyr and 
NR Lyr, taking the standard deviation as a measure of the 
uncertainty. In this way we arrived at the following Kp-V 
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Figure 8. Four panels comparing Fourier parameters (V-band) for the Kepler non-Blazliko ab-type RR Lyr stars (large black squares) 
and for 24905 RR Lyr stars in the central regions of the Large Magellanic Cloud. The LMC data are from the Soszynski et al. (2009) 
OGLE-III study, with colour and symbol coding as follows: red filled dots (17693 RRab stars); blue filled squares (4957 RRc stars); green 
open triangles (986 RRd stars), and purple open diamonds (1269 RRe stars). Because of the limited x- and y-ranges not all of the LMC 
RR Lyr stars are represented in the graphs. 



transformations: 



Atot{V) = Atat{Kp) + (0.14 ± 0.01), 

RT{V) = RT{Kp) + (0.002 ± 0.006), 

Ai{V) = Ai{Kp) + (0.054 ± 0.002), 

AsiV) = AsiKp) + (0.018 ± 0.003), 

R21 (Y) = R21 {Kp) - (0.013 ± 0.004) , (2) 

RsiiV) = Rzi{Kp) ~ (0.004 ± 0.009), 

(t>2i{V) = <t>2iiKp) - (0.089 ± 0.021), 

<l>3i{V) = 031 (A» - (0.151 ±0.026). 



The most uncertain of these equations are those for RT and 
i?3i . While the percent uncertainties are large the offsets are 
small; and besides, these two quantities are rarely used in 
correlations to derive physical properties. 

A comparison of the ASAS Atot{V) values (given in Ta- 
ble 2, column 7) with the Kp values (given in column 8 of 
Table 1) is not helpful since several of the stars show larger 
than expected differences (for example, V838 Cyg). The ac- 
curacy of the above transformation equations is expected 
to improve with future BVI photometry. In §6 these offset 
equations are used to convert the well-established T/-band 



correlations relating Fourier parameters and physical char- 
acteristics to Kp correlations. 

5.3 Kepler vs. globular cluster RR Lyr Stars 

In Figure 7 the Fourier parameters for the 19 non-Blazhko 
stars in the Kepler field (derived using the Kp photometry 
but transformed to V values using the offsets given in Eq. 2) 
are compared with parameters derived from V^-photometry 
for 177 RR Lyr stars in several well-studied globular clusters 
(from Kovacs & Walker 2001). The same symbols and colour 
coding are used in all four panels. 

The upper-left panel is analogous to the P-Atot dia- 
gram and can be used to derive metallicities (see §6.1 be- 
low). The two diagonal lines represent mean relations when 
the globular cluster data (small black dots) are sorted into 
two [Fe/H]BB bins (here 'BB' refers to the Butler-Blanco sys- 
tem, which approximates the Carretta-Gratton system - see 
below): a metal-poor bin consisting of 19 stars (surrounded 
by blue squares) with metallicities between —1.70 and —1.99 
dex and average —1.8 dex (BB-scale); and an intermediate- 
metallicity bin consisting of 39 stars (circled in red) with 
[Fe/H]BB between —0.97 and —1.23 dex and average —1.1 
dex. The equations of the lines are: = 5.556 logP-f 0.2920 
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(upper metal-poor bin) and (l>?n ~ 6.200 logP + 3.615 
(lower intermediate metallicity bin). Four of the ifep/er non- 
Blazhko stars are clearly richer than [Fe/H]BB=— 1.1 dex, 
while the remainder are apparently more metal poor. 

The axes of the other three panels in Fig. 7 are similar 
to those in Fig. 6. The upper-right panel shows that the Ke- 
pler stars appear to be drawn from a distribution similar to 
that of the globular clusters. In particular, the Kepler stars 
with low R31 values are not unusual, except possibly that 
they all have relatively high R21 values. Since the globular 
cluster stars of higher metallicity (red open circles) all tend 
to reside on the right side of the diagram this separation 
is probably a metallicity effect, supporting our conclusion 
that V784 Cyg, V782 Cyg, KIC 6100702 and V2470 Cyg 
are metal rich. Likewise, the lower-left panel shows that the 
majority of the Kepler stars do not differ from the stars in 
globular clusters. Note too that there is very little metallic- 
ity discrimination in this plane. The two stars located at the 
extreme upper edge of the envelope of the globular cluster 
distribution are AW Dra and V784 Cyg. Finally, the lower- 
right panel shows close agreement between the phase param- 
eters of the Kepler and globular cluster RR Lyr stars, which 
supports the conclusion drawn earlier (lower-right panel in 
Figure 6) that there is a strong approximately-linear corre- 
lation between 021 and ^31. The diagram also shows very 
little dependence on metallicity. 

5.4 Kepler vs. LMC field RR Lyr stars 

In Figure 8 we compare the R21, R31, <j)2i said 4)^i Fourier 
parameters for the 19 ifepZer non-Blazhko RR Lyr stars with 
those for the field RR Lyr stars in the central regions of the 
Large Magellanic Cloud. The LMC data are from the mas- 
sive OGLE-III Catalogue of Variable Stars by Soszynski et 
al. (2009), which comprises almost 25000 RR Lyr stars. All 
the parameters are y-band values, the LMC values hav- 
ing been transformed from I to V using the transformation 
equations given by Morgan, Simet & Bargenquast (1998), 
and the values for the Kepler RR Lyr stars from Kp to V 
using the Eq. 2 offsets. 

The upper-left panel shows that the Kepler stars have 
R21 values in the narrow range 0.43 to 0.51, which is much 
smaller than the range seen for the LMC stars, 0.25 to 0.55. 
In particular, none of the Kepler RRab stars are among the 
LMC RRab stars in the '_R2i dropdown' at logP ~ —0.19, 
I.e., P ~ 0.65 d. Note that the upper-left panel of Fig. 6 can 
be used to identify individual Kepler stars. 

The R21 and i?3i amplitude ratios are compared in the 
upper-right panel of Fig. 8. This graph, which is similar 
to the upper-right panels of Figs. 6 and 7 (and again. Fig. 6 
can be used to identify individual if epier stars), further puts 
into perspective the small ranges of 7?2i and Psi seen for the 
if^ep/er RRab stars, both relative to that for the LMC RRab 
stars, and relative to the ranges seen for the other types of 
RR Lyr stars (RRc, RRd and RRe). 

The bottom panels of Fig. 8 are both 'metallicity di- 
agnostic' diagrams, such as those discussed in §6.1 below. 
We have already seen (upper-left panel of Fig. 7) that the 
Kepler stars sort into four metal-rich and 15 metal-poor 
stars, that the metal-poor stars are similar to the RRab 
stars found in both Oosterhoff type I and type II galactic 
globular clusters, and that the metal-rich stars have no coun- 



terparts among the Kovacs & Walker (2001) sample of glob- 
ular cluster RR Lyr stars. Here we see that both the metal- 
rich Kepler stars (the clump of four stars at logP~— 0.3 and 
02i~4.3) and the metal-poor ifep/er stars have counterparts 
in the LMC, and that the metal-rich stars appear to have 
relatively long periods for metal-rich RRab stars (i.e., the 
subgroup of LMC RRab stars with shorter periods at a given 

031 ). 



6 PHYSICAL CHARACTERISTICS 

RR Lyrae stars provide fundamental insight into the late 
evolution of low mass stars, and in particular the instabil- 
ity strip (IS) region of the horizontal branch (HB). Being 
present in many globular clusters (and galaxies) exhibiting 
different HB typefl their properties are related to the evo- 
lutionary and chemical histories of these systems. Thus it is 
desireable to derive (from photometry, spectroscopy and de- 
tailed pulsation and stellar evolution models) their physical 
characteristics, such as the mass M, luminosity L, effective 
temperature Tetr, iron abundance [Fe/H], chemical composi- 
tion (X,Y,Z), age, etc. An equally important goal is to under- 
stand the interdependencies of these quantities. Significant 
early contributions are summarized in the papers by Christy 
(1966), Iben (1971), Smith (1995), Sandage & Tammann 
(2006, hereafter ST6), and Catelan (2009). This knowledge 
becomes critical when RR Lyr stars (and Cepheids) are used 
to derive accurate distances within and beyond our Galaxy 
(see Tammann, Sandage & Reindl 2008, hereafter TSR8; 
Sandage & Tammann 2008, hereafter ST8). 

Of course high-dispersion spectroscopy permits direct 
measurements to be made of some of these quantities; how- 
ever, such observations are time-consuming and often im- 
practical. Thus efforts have been on-going to derive physi- 
cal characteristics from photometry only, in particular from 
pulsation periods, mean magnitudes and colours, quanti- 
ties that characterize the shapes of the light-curves {e.g., 
amplitudes, risetimes, Fourier phase parameters), and most 
recently, from fitting detailed light curve shapes to non- 
linear convective pulsation models (see Bono, Castellani 
& Marconi 2000, Marconi & Clementini 2005, Marconi & 
Degl'Innocenti, 2007). 

A potential difficulty for the present Fourier investi- 
gation is that almost all the available calibration relations 
use empirical correlations derived from V-band photometry, 
while for the stars studied here we have Kp photometry, and 
only a limited amount of high-precision V photometry and 
colour information. Furthermore, the mean Kp magnitudes 
from the KIC are somewhat uncertain for giant stars. De- 
spite these limitations we have shown in the previous section, 
at least for AW Dra, FN Lyr and NR Lyr, three RR Lyr stars 
with quite different periods and light curve shapes, that dif- 
ferences between the V and Kp Fourier parameters appear 
to be small and systematic. For this reason we have chosen 
to present approximate physical characteristics based on the 

^ The HB type of a globular cluster is defined as 
(B-R)/(B+V+R), where B, R and V are the number of 
blue HB stars in the cluster, the number of red HB stars, and 
the number of RR Lyr stars, respectively (Lee 1989). A table of 
values is given in Lee, Dcmarque &: Zinn (1994). 
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extant IZ-band correlations applied to Kp correlations trans- 
formed using the observed Kp-V parameter offsets. Our re- 
sults are summarized in Tables 4-6 and their derivation is 
discussed next. 



6.1 Iron Abundances 

6.1.1 Background: [Fe/H] from P-shifts in P-A diagram 

Starting with the papers by Oosterhoff (1939, 1944), Arp 
(1955) and Preston (1959), the period-amplitude {P-A) di- 
agram, also known as the 'Bailey diagram', has been found 
to be a useful tool for deriving iron-to-hydrogen ratios for 
field and cluster/galaxy RR Lyr stars. The amplitude usu- 
ally employed is the B ot V total amplitude, ^tot, and often 
the abscissa is logP. The basic P-j4tot-[Fe/H] correlation is 
such that for a given value of Atot lower- met allicity stars 
tend to have longer periods than higher-metallicity stars. 
Sandage (2004, hereafter S04) has called this the 'OAP 
period-metallicity correlation'. Other recent discussions of 
the P-A diagram and period shifts are given by Di Cri- 
scienzo et al. (2004), Bono et al. (2007, hereafter BCD7) 
and Sandage (2006 and 2010, hereafter S06 and SIO). 

The basic explanation for P-shifts in the Bailey dia- 
gram follows from the 'stacked HB luminosity levels' model 
first proposed by Sandage (1958), the latest refinement of 
which is given in SIO (see his figs. 1 and 3). Owing to the 
Ritter (1879) pulsation relation, Py/p — Q (where Q is a 
constant over a wide range of L, Tcfi and [Fe/H]), lines of 
constant density (p), and therefore period, cut diagonally 
across HBs of different luminosity (see fig. 13 of Sandage, 
Katem & Sandage 1981, hereafter SKS). When coupled with 
a monotonic Atot-Tes correlation (in the sense that Atot in- 
creases with distance from the red edge of the IS - see, for 
example, fig.l2 of SKS, and fig. 12 of Di Criscienzo et al. 
2004) one expects to see at a given Atot more luminous HB 
stars having longer periods. In fact, P-shifts and correlations 
such as the logP-Atot-[Fe/H] correlation are expected in all 
diagrams where the abscissa is pulsation period and the or- 
dinate is a light curve descriptor that correlates with Atot. 
Two such descriptors are the risetime (RT) and the Fourier 
4>3i phase parameter (defined in §4.2) - see Fig. 10 below. 
For example, in the RT vs. logP diagram the P-shift at a 
given RT between M3 and M15 amounts to AlogP = 0.055 
(see fig. 10 of SKS). Period shifting can also be seen between 
M3 and M2 in the 031 vs. period diagram (see fig. 2 of Ju- 
rcsik et al. 2003), and between metal-rich and metal-poor 
field RR Lyr stars (see fig. 3 of S04). Thus, such correlations 
have the potential for measuring [Fe/H] (S04; Kovacs 2005; 
Cacciari, Corwin & Carney 2005, hereafter CCC5; S06,S10; 
Sodor, Jurcsik & Szeidl 2009). 

Identifying the less luminous HB in the Sandage model 
with that of the unevolved zero-age horizontal branch 
(ZAHB) stars in the IS of a relatively metal-rich system, 
and the more luminous HB with the same near-ZAHB stars 
in the IS of a more metal-poor system, leads directly to 
the OAP period-metallicity correlation. Indeed, ZAHB stel- 
lar evolution models clearly show that the mean luminosity 
levels in the IS for metal-poor ZAHB stars are higher than 
those for more metal-rich ZAHB stars {e.g., Faulkner 1966; 
Faulkner & Iben 1966; Iben & Faulkner 1968; Hartwick, 
Harm & Schwarzschild 1968; Iben & Rood 1970; Iben 1971; 



Sweigart & Gross 1976; Lee, Demarque & Zinn 1990,1994, 
hereafter LDZ90 and LDZ94; Dorman 1992; VandenBerg et 
al. 2000; Demarque et al. 2000, hereafter DOO; Cassisi et al. 
2004; Pietrinferni et al. 2004, 2009; Dotter et al. 2007, etc.). 
Indeed, in the IS at logTcH = 3.85 the M15 HB is ~0.22 
mag brighter than that for M3 (SKS). And, for [Fe/H] val- 
ues ranging from —2.3 to —0.5 the mean absolute magni- 
tudes for the RR Lyr stars range from {Mv(RR)) ~ 0.40 to 
0.85 mag, with, for all but the most extreme clusters, little 
dependence on HB type (see figs. 1,2 of DOO; and fig. 3 of Cas- 
sisi et al. 2004). The models also show that in the instability 
strip at the luminosity level of RR Lyr stars the variation 
in L(ZAHB) is due primarily to a difference in mean mass, 
with (Ai) being higher in more metal-poor system^. 

6.1.2 Period shifts caused by post-ZAHB evolution 

Despite the general success of the Bailey diagram as a tool 
for deriving [Fe/H] values it has become apparent since the 
early 1980's that simple correlations of period-shift with 
metallicity do not exist in all cases and factors other than 
metal abundance play a role in the interpretation of these 
diagrams and correlations. Chief among these effects is evo- 
lution away from the ZAHB. 

All HB stellar evolution models show tracks leading 
away from the ZAHB toward the asymptotic giant branch, 
with HB lifetimes ~50-100 Myr (core helium burning phase). 
While the evolution can be either redward or blueward it 
almost always is in the direction of higher luminosity, with 
the pace quickening as the evolution progresses. As L in- 
creases the period becomes longer and larger P-shifts are 
seen in the P-A diagram. However, most HB stars are ex- 
pected to be located near to the ZAHB, with fewer stars 
in advanced evolutionary phases (CCC5 put the fraction of 
evolved RR Lyr stars in M3 at ~14%). Indeed skewed HB 
and RR Lyrae luminosity distributions are observed in many 
globular clusters {e.g., fig. 15 of Sandage 1990, hereafter S90). 
Such skewed distributions are also seen in simulations of HB 
and RR Lyr stars (see fig. 10 of Marconi et al. 2003, and fig. 3 
of Catelan 2004), and in period shift diagrams (see fig. 6 of 
Kunder & Chaboyer 2009). 

As a first step towards identifying and correcting 
for such evolutionary effects Sandage (1981a,b; hereafter 
S81a,b) introduced the concept of 'reduced' period, logP' 
= logP-f0.336A7nboi, where Amboi = rn^oi minus the 
mean apparent magnitude of the bulk of the RR Lyr stars 
in a given cluster, and the multiplicative factor follows 
from the van Albada & Baker (1971,1973) version of the 
Pyfp = Q relation, and demonstrated the existence of a 
period-luminosity-amplitude relation for equal metallicity 
RR Lyrae stars. The net effect of making this correction 
is to reduce the period shifts of the evolved RR Lyr stars, 
resulting in tighter correlations in the Atot vs. logP' dia- 
gram. In this way S81a,b identified several stars in Messier 3 

^ Based on the Atot-ToS correlation one expects also to see pe- 
riod shifting in the logTgff vs. logP diagram, and indeed fig. 11 
of SKS shows at 'any given temperature' an observed shift be- 
tween M3 and M15 of AlogP=0.070. However, when comparing 
model predictions with observations one must 'let the tempera- 
ture be cooler for more metal-poor RR Lyr stars at their higher 
luminosities' (Simon & Clement 1993, S93, S06). 
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that are extra-luminous with longer periods at a given Ab 
than the bulk of the RRab starfl. Many other globular 
clusters are now known to contain extra-luminous large- 
period-shift stars, including 47 Tuc (V9 - see Storm et al. 
1994), NGC 6388 and NGC 6441 (Pritzl et al. 2000, 2001, 
2002; Catelan et al. 2006), Messier 5 (Kaluzny et al. 2000, 
CCC5), NGC 5896 (Alves et al. 2001), IC 4499 (V9 and V54 
- see CCC5; Kunder et al. 2011; Walker et al. 2011), and 
Messier 13 (Sandquist et al. 2011). 

The reduced-period concept is useful when deriving 
mean metallicities for ensemble populations of RR Lyr stars 
in globular clusters and galaxies {e.g., Dall'Ora et al. 2003), 
and for field stars that are near the ZAHB (which are as- 
sumed to be in the majority), but the method fails for 
evolved RR Lyrae stars: what looks in the Bailey diagram 
like a metal-poor ZAHB star might well be an evolved metal- 
rich star, i.e., an evolved Ool star masquerading as an un- 
evolved OoII star. 

Thus, the amount of P-shifting is determined by the 
amount of L-evolution, which in turn depends on the the HB 
morphology of the parent population (LDZ90, Lee 1990j3- 
When RR Lyr stars originate from a position on the ZAHB 
that is within the IS, which occurs in Ool clusters with HB 
type near zero, such as M3 and M4 with HB types (from 
LDZ94) of 0.08 ± 0.04 and -0.07 ± 0.10, respectively, the 



^ Sandagc's candidate evolved stars in M3 include V24, V60, V65, 
V96 and V124, all of which can be seen (labelled) in the S81a fig.2 
"^bol VS. logTeff diagram and in the S81a fig. 3 period-amplitude 
diagrams. Subsequent studies of MS by Kaluzny et al. (1998) and 
Clement &; Shelton (1999) confirmed V65 as an evolved star and 
also have identified V14 and V104 as post-ZAHB stars. Post- 
ZAHB evolution in MS was also investigated by Marconi et al. 
(2003) and by Jurcsik et al. (2003). Jurcsik et al. used "luminosity 
bins" and Fourier parameters and found that RR Lyr "stars close 
to the ZAHB show Ool type properties, while the brightest stars 
have Oo II statistics regarding their mean periods and RRab/RRc 
number ratios" , confirming the earlier result of Clement & Shel- 
ton (1999). Thus, again, the Gosterhoff dichotomy seems to be 
connected with evolutionary effects. More recently, CCC5 con- 
firmed V24 and V14 as evolved stars and identified three more 
post-ZAHB stars based on their long periods and relatively high 
luminosities (VS, VS5, V67). Most recently, Valcarce & Cate- 
lan (2008) identified in MS ~15 non-Blazhko and another ~10 
Blazhko RR Lyr stars in an advanced post-ZAHB evolutionary 
state (see their Fig.2). 

^ The main parameter controlling the distribution of stars along 
the HBs of globular clusters (GCs) is metal abundance (Sandage 
& Wallerstein 1960), with metal-rich clusters (e.g., NGC 6S56, 
47 Tuc) having predominantly red HBs and metal-poor systems 
(e.g., M15, M9) having predominantly blue HBs. Since the discov- 
ery by Sandage & Wildey (1967) and van den Bergh (1967) that 
there is not a unique correlation between [Fe/H] and HB mor- 
phology, and the recognition that the GCs with the bluest HBs 
are not the most metal poor clusters (Renzini 1983), the search 
has been on for the so-called 'second parameter'. The two leading 
candidates are differences in age (see Rood &; Iben 1968; Lee, De- 
marque & Zinn 1994; Lee et al. 1999, 2007; Rey et al. 2001; Dotter 
et al. 2010), and variations in the initial helium abundance (see 
Sandage & Wildey 1967, van den Bergh 1967, Hartwick 1968, 
Busso et al. 2007). Other candidates include variations in the 
amount of mass loss along the red giant branch (RGB), rotation 
rates, and the amount of helium mixing deep inside progenitor 
RGB stars (see Sweigart & Catelan 1998). 



evolution manifests itself in the IS region of the HR dia- 
gram as a vertical widening of the HB with RR Lyr stars of 
similar mass spread over a range of L and TaS (see S90; and 
fig.l of Jurscik et al. 2003, hereafter J03; and fig.2 of Val- 
carce & Catelan 2008) . On the other hand, when the original 
location on the ZAHB was blueward and outside of the IS, 
which occurs in OoII clusters with HB types that are large 
and positive, such as M15 (see Sobeck et al. 2011), M92 
(see Roederer & Sneden 2011) and M2 (Lee & Carney 1999) 
with HB types of0.72±0.10, 0.88±0.08 and 0.96±0.10 (all 
from LDZ94), respectively, most if not all the RR Lyr stars 
are evolved and few if any lower-luminosity ZAHB stars are 
present (see LDZ90, DOO). In the case of M2 Lee & Car- 
ney (1999) concluded that its RR Lyrae stars are ~0.2 mag 
more luminous than those in M3, owing to all the M2 stars 
having evolved away from the blue side of the HB, whilst 
the bulk of the M3 stars lie near the ZAHB. Additional ob- 
servational support for this idea comes from the clear cor- 
relation of metallicity and HB intrinsic vertical width, with 
more metal-rich systems having greater vertical widths (see 
Fig. 16 of S90). Schematic tracks (along with grid lines of 
constant period and amplitude) comparing the post-ZAHB 
evolution of stars in M2 (and M15) with the tracks for M3 
(and NGC6441) are given in Fig. 3 of SIO. Completing the 
picture, in clusters with red HBs {i.e., HB types that are 
large and negative, such as 47 Tuc and M107, with HB types 
of —1.0 ± 0.03 and —0.76 ± 0.08, respectively) the evolution 
occurs mainly redward of the IS and few, if any, RR Lyr 
stars are expected (exceptions to this rule include systems 
with red HBs and extended blue HBs, such as NGC 6388 
and NGC 6441, - see references above). Thus, the luminos- 
ity of an individual RR Lyr star in a given cluster depends 
not only on its metallicity (usually assumed to be that of 
the parent cluster) but on its evolutionary state, which in 
turn is determined by the HB-type of the parent cluster and 
the amount of evolution away from the ZAHB. 

In closing, a potential solution to the evolution degener- 
acy problem in the Bailey diagram is to identify the evolved 
RR Lyrae stars using some sort of light curve index, such 
as the JK96 'compatibility criterion' (see Kovacs & Kanbur 
1998, Clement & Shehon 1999, Nemec 2004, CCC5). Of par- 
ticular interest in this regard is fig. 15 of CCC5, which shows 
large P-shifts for several stars in the Ab vs. \ogP diagram, 
but no such offset in the ^31 vs. \ogP diagram; whether post- 
ZAHB stars stand out more in the P-A diagram than in the 
period-(?!>3i diagram remains an open question. In any case it 
is important to excercise caution when using the individual 
Fourier-based metallicities derived here. 

6.1.3 Metallicity Scales 

For the Kepler non-Blazhko RRab stars we have used the 
correlations involving all three light curve descriptors to es- 
timate [Fe/H] values (see §6.1.4). The resulting metallicities 
are reported in Table 4. The metal abundances based on 
the S04 equations are on the Zinn & West (1984, hereafter 
ZW) scale, and those based on J98's equation 1 are on 
the Carretta & Gratton (1997, hereafter CG) scale. For 
comparison purposes the latter (i.e., the cjili-metallicities) 
have been transformed to the ZW scale and to the Carretta 
et al. (2009, hereafter C9) scale. Although several other 
metallicity systems have been proposed {e.g., Layden 1994, 
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Figure 9. Period-amplitude and three other metallicity diagnostic diagrams for the Kepler non-Blazhko RRab stars. For all four graphs 
the abscissa is logP, all the stars have been labelled, and the four candidate metal-rich stars are plotted as open squares. AW Dra, 
FN Lyr and NR Lyr arc plotted twice in each panel, with the points derived from the Kp photometry (solid black squares) connected 
to the points derived from the high-precision y-band photometry (open diamonds) - the consistency of the Kp-V offsets seen here (and 
in Figs. 6 and 10) is quite striking. The diagonal dashed lines, representing the relations for the M3 RR Lyr stars (assumed to have 
[Fe/H]zw = —1.6 dex), are from Equation 5 (S04eq6, top left). Equation 6 (S04eq7, bottom left) and fig. 15 of CCC5 (bottom right). 



Jurcsik 1995, Butler & Blanco - see S04, and Kraft & Ivans 
2003) , the C9 scale is most extensive. It is based on a homo- 
geneous data set derived from ~2000 modern spectra, and is 
gaining considerable favour - see, for example, its adoption 
by Harris (1996) for his on-line globular cluster catalog at 
http://www.physics.mcmaster.ca/^harris/mwgc.dat |. In 
addition, equations are provided for transformation between 
the various scales. When transforming between the ZW and 
CG systems the following equation (from S04, footnote 1) 
has been used: [Fe/H]zw = 1-05 [Fe/H]cG - 0.20. And when 
transforming from the ZW and CG scales to the C9 scale 
we have used (from C9): [Fe/H]c9 = 1-137 [Fe/H]cG - 0.003; 
and [Fe/H]c9 = -4.13 + 0.130 [Fe/H]zw - 0.356 [Fe/H]|w. 



6.1.4 Metallicities for the Kepler stars 

Figure 9 contains, for the Kepler non-Blazhko RRab stars, 
four such 'metallicity diagnostic diagrams'. The principal 
observation here is the similarity, to first order, of the dis- 
tribution of the points in the four graphs. The top-left panel 
is the classical P-Atot diagram, with total Kp amplitude 
plotted along the ordinate. Traditionally, lines of constant 
[Fe/H] serve as metallicity calibration lines; these run diag- 
onally from the upper-left to the lower-right corner (as seen 
in three of the panels). 

The top-right panel has A\ along the ordinate. As 
expected, because Ai and Atot are highly correlated (see 
Fig. 2), the distribution of the points closely resembles that 
in the Atot vs. logP diagram. In principle the small range 
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of Ai (only ~0.2 mag) compared with that of Atot (nearly 
0.7 mag here, but often larger) makes the P-Ai diagram less 
useful for estimating [Fe/H]; but with such high-precision 
data as provided by Kepler the derived [Fe/H] results should, 
in principle, be similar. 

The bottom panels of Fig. 9 show risetime (left) and 
031 (right) plotted against logP, where both RT (reverse 
scale) and (reverse scale) increase with increasing logP. 
The bottom-right panel is the ^Kepler non-Blazhko RRab 
only' version of the top-left panel of Fig. 7, and also is sim- 
ilar to Fig. 11. The dashed line in the bottom-right panel 
(from fig. 15 of CCC5), the equation of which is (jj^i — 
3.124 -|- 5.1281ogP, represents the observed V-band relation 
for the RRab stars in M3 - it is encouraging that the slope of 
the line appears to agree with the slope of the Kepler metal- 
poor stars. According to this calibration NR Lyr is the most 
metal-poor star in the sample, followed by KIC 7021124 and 
FN Lyr, and most of the RRab stars have metallicities simi- 
lar to those in MS, which has [Fe/H]= -1.50±0.05 dex (Car- 
retta et al. 2009 scale; W. Harris' on-line catalog). Metallic- 
ities for the other stars, including the four stars suspected 
of having high [Fe/H] values, are discussed below. 

Figure 10 contains the RT vs. Atot diagram (left panel) 
and the (f>3i (which is vr larger than (p^i) vs. Atot diagram 
(right panel) for the Kepler stars. In general, as the total 
amplitude increases, both RT and (pf^i decrease, with linear 
fits described by the equations: RT = — O.llSyltot + 0.256 
and (f)fji = — 1.75yltot + 6.75; however, owing to the large 
scatter about these lines, and the suspected large range 
in metallicities of the Kepler stars, we believe these equa- 
tions to be unreliable (and for this reason they were not 
plotted). Theoretical models (see the middle-left panels of 
Figs. 14 and 15 below) suggest that it is more probable that 
there are families of lines (or curves), one for each metal- 
licity. Unfortunately, the present observational data are in- 
sufficient to establish these. It is noteworthy that AW Dra, 
and to a lesser extent V784 Cyg and V894 Cyg, have larger 
(j}fji values and larger RTs than other stars of comparable 
^tot. And, while it appears that all the low-amplitude stars 
have large RTs and large (l>3i values, this conclusion is based 
on only four low-amplitude candidate metal-rich stars and 
two low-amplitude long-period candidate metal-poor stars 
(KIC 7030715, AW Dra), one of which (AW Dra) is a possi- 
ble outlier in both graphs and which may be a long-period 
evolved star of intermediate-metallicity (if our preliminary 
CFHT spectra result is correct - see §6.1.5 below) rather 
than a near-ZAHB star of low-metallicity. In any case, both 
KIC 7030715 and AW Dra have considerably longer periods 
than the four high-metallicity stars. 

In Figure 11 the Kepler data are compared with the 
galactic field star data analyzed by S04 (see his fig. 3). The 
field RR Lyr stars are the same as those studied by Simon 
& Lee (1981), Simon & Teays (1982) and Simon (1988), the 
details of which are given in table 1 of S04. Fig. 11 is com- 
parable to the bottom-right panel of our Fig. 9, which shows 
only the /fep/er non-Blazhko stars. Because the Simon et al. 
stars have known [Fe/H]zw values from Layden (1994) it is 
possible to sort them by metallicity. Following Sandage we 
plotted two subsets of the data, which have average metal- 
licities ([Fe/H]) = -0.35 and -1.63 dex (ZW scale); the 
diagonal lines were obtained by substituting these values 
into eq. 3 of S04. This diagram is also the galactic field star 



version of the upper-left panel of Figure 7, which included 
RR Lyr stars in galactic and Magellanic Cloud globular clus- 
ters. The chief difference is that the calibration has now been 
extended to more metal-rich stars. 

Fig. 11 shows that the stars KIC 6100702, V2470 Cyg, 
V782 Cyg and V785 Cyg, all of which have long risetimes, 
are solidly in the metal-rich region of the diagram, with 
the highest (j>3i values and all with periods shorter than 
0.55 d. The other Kepler stars are a better match to the 
field stars that have [Fe/H]zw=— 1.63 dex. AW Dra and 
KIC 7030715, with long periods, intermediate risetimes, and 
high (^31 values, are seen to have metallicities ~ —1.64 dex 
(ZW scale). Three of the four metal-rich stars (V782 Cyg, 
KIC 6100702 and V2470 Cyg) have [Fe/H]zw = -0.4 dex, 
and V784 Cyg is the most metal-rich star in our sample with 
[Fe/H]zw = —0.2 dex. The most metal-poor stars are seen 
to be more metal-poor than [Fe/H]zw = —1.6 dex. 

Equation 1 of J98, which is identical to equation 3 of Ju- 
rcsik & Kovacs 1996, describes lines of constant [Fe/H]cG in 
the period-03i(l/) diagram. Using the offset formula for <^3i 
given in Eq. 2 the Kp version of the J98 equation becomes 

[Fe/H]cG = -5.241 - 5.394 P-f 1.345 <?!>ii(/Cp). (3) 

Metallicities calculated with this formula are listed in col- 
umn 2 of Table 4 under 'J98eql'. For comparison purposes 
the corresponding values on the ZW and C9 scales also are 
listed. 

Our second estimate of [Fe/H] is derived from the 
[Fe/H]zw - logP - (^31 relation of S04 (his eq. 3 - note that 
the cosine version of the Fourier phase parameter is em- 
ployed here) - see lower-right panel of Fig. 9. Applying the 
(^31 offset given in Eq. 2, Sandage's relation becomes 

[Fe/H]zw = -6.217- 7.012 logP-f 1.411 <j!.;ii(A'p), (4) 

where the intercept term has a total uncertainty of ±0.023, 
and the (p^i and logP coefficients have respective uncertain- 
ties of ±0.014 and ±0.071 (from S04). Metallicities derived 
from Eq. 4 are listed in Table 4 under 'S4eq3'. 

Implementation of the P-A diagram (see upper-left 
panel of Fig. 9) for deriving metallicities was via eq. 6 of 
S04, combined with the Kp-V offset for Aot (Eq. 2): 

[Fe/H]zw = -1.453 AtotiKp) - 7.990 logP - 2.348, (5) 

where the errors in the coefficients (from S04) are ±0.027, 
±0.091 and ±0.043, respectively. The derived [Fe/H]zw val- 
ues are listed in Table 4 under 'S4eq6'. According to this 
formula the most metal-poor star is AW Dra, 0.2 dex more 
metal-poor than NR Lyr, and 0.6 dex more metal-poor than 
the estimates derived form the (^31 equations! Given its lo- 
cation in the (j>^i vs logP diagram, and our recently derived 
spectroscopic metallicity of — 1.33±0.08 dex, it is our opinion 
that this is an evolved RR Lyrae star; curiously, it follows 
the pattern of evolved RR Lyr stars seen in M3 by CCC5 
(their fig. 16). The metal-rich stars also are predicted to be 
more-metal poor than the values given by the (j>3i estimates. 

Our fourth set of [Fe/H] estimates was based on the RT 
vs. logP diagram (lower- left panel of Fig. 9). The metallic- 
ities were calculated using Sandage's relation between RT, 
P and [Fe/H] (his eq. 7), combined with our V-Kp offset 
(Eq.2): 

[Fe/H]zw = 6.33 KT{Kp) - 9.11 logP - 4.59, (6) 
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Figure 10. Risetime vs. total amplitude diagram (left panel) and Fourier phase parameter (f)^^ vs. total amplitude diagram (right panel) 
for the Kepler non-Blazhko RRab stars. The symbols are the same as in the bottom-right panel of Fig. 9 (i.e., the four metal-rich stars 
have been plotted with open squares), and again AW Dra, FN Lyr and NR Lyr are plotted twice in each panel, indicating the magnitudes 
and directions of the (systematic) V-Kp offsets. 



Table 4. Metal abundances based on the Kp photometry and the 
equations of Jurcsik (1998) and Sandage (2004). 



Star 












[Fc/H] 




















^tot 


RT 






J98eql 






S4eq3 


S4eq6 


S4eq7 






CO 


, ZW, C9 




ZW 


ZW 


ZW 


(1) 






(2) 






(3) 


(4) 


(5) 


NR Lyr 


-2 


04, 


-2.34, - 


-2 


32 


-2.28 


-2.13 


-2.16 


V715 Cyg 


-1 


18, 


-1.44, - 


-1 


34 


-1.45 


-1.17 


-0.74 


V782 Cyg 


-0 


25, 


-0.47, - 


-0 


29 


-0.51 


-0.86 


-0.76 


V784 Cyg 


+0 


06, 


-0.14,- 


fO 


07 


-0.18 


-1.09 


-0.87 


KIC 6100702 


-0 


14, 


-0.35, - 


-0 


17 


-0.38 


-0.69 


-0.48 


NQ Lyr 


-1 


55, 


-1.83, - 


-1 


77 


-1.86 


-1.68 


-1.57 


FN Lyr 


-1 


62, 


-1.90, - 


-1 


84 


-1.94 


-1.70 


-1.31 


KIC 7021124 


-1 


79, 


-2.08, - 


-2 


04 


-2.08 


-1.91 


-1.83 


KIC 7030715 


-1 


39, 


-1.66, - 


-1 


58 


-1.60 


-1.97 


-1.96 


V349 Lyr 


-1 


46, 


-1.73, - 


-1 


66 


-1.77 


-1.47 


-1.08 


V368 Lyr 


-1 


27, 


-1.53, - 


-1 


44 


-1.53 


-1.27 


-0.64 


V1510 Cyg 


-1 


56, 


-1.83, - 


-1 


77 


-1.86 


-1.89 


-1.57 


V346 Lyr 


-1 


55, 


-1.83, - 


-1 


76 


-1.86 


-1.84 


-1.52 


V350 Lyr 


-1 


56, 


-1.84, - 


-1 


78 


-1.86 


-1.96 


-1.62 


V894 Cyg 


-1 


51, 


-1.79, - 


-1 


72 


-1.82 


-2.01 


-1.46 


V2470 Cyg 


-0 


48, 


-0.71, - 


-0 


55 


-0.74 


-1.13 


-0.99 


V1107 Cyg 


-1 


30, 


-1.56, - 


-1 


48 


-1.60 


-1.56 


-1.38 


V838 Cyg 


-1 


30, 


-1.56, - 


-1 


48 


-1.58 


-1.40 


-0.84 


AW Dra 


-1 


47, 


-1.74, - 


-1 


67 


-1.68 


-2.34 


-2.06 



where the S04 intercept has been increased by 0.01 (repre- 
senting the RT conversion from V to Kp). Application of 
this equation gave the metallicities listed in Table 4 under 
'S4eq7'. 



6.1.5 Discussion of the derived metallicities 

In the absence of spectroscopic metal abundance^ which 
of the [Fe/H] values in Table 4 are the most reliable? One 
area of potential bias comes from the fact that all the metal- 
licity correlations used in §6.1.4 assume that the stars are 
unevolved. This probably is true for most of the Kepler non- 
Blazhko stars; however, there may a few evolved stars in the 
sample, the most likely candidate being AW Dra. In such a 
case the period will have been increased by luminosity evo- 
lution and the derived [Fe/H] will be systematically smaller 
than the true value. Another potential bias occurs for the 

* High-resolution echelle spectra of the brightest 18 Kepler 
RR Lyr stars (non-Blazhko and Blazhko, RRab and RRc) recently 
have been acquired with the ESPaDOnS spectrograph on the 
Canada-France-Hawaii (CFHT) 3.6-m telescope (Nemec, Ripcpi, 
Chadid et al. 2011, in preparation). These spectra currently are 
being measured and preliminary informative results are available 
for two of the non-Blazhko RRab stars: (1) For KIC 6100702 a 
spectroscopic metal abundance of —0.18 it 0.06 dex has been do- 
rived - this compares very favourably with the two 03i-based 
metallicities, [Fe/H]c9 = -0.17 dex and [Fe/H]zw = -0.38 dex 
(see columns 2 and 3 of Table 4), both of which are more metal- 
rich than the Atot and RT values. Thus for one of the four sus- 
pected metal-rich non-Blazhko RRab stars there is spectroscopic 
confirmation of high metallicity; and (2) for AW Dra a spectro- 
scopic [Fc/H]= —1.33 ±0.08 dex has been derived. If correct this 
is more metal-rich than any of the derived photometric estimates. 
Our suspicion is that AW Dra, which exhibits a large period-shift 
in the Atot vs. logP diagram and in the Ai vs. logP diagram (see 
upper two panels of Fig. 9), and is an outlier in the vs. Ai 
diagram and in the </>|]^ vs. (p^i diagram (see bottom two panels 
of Fig. 6), and an outlier in both Fig. 10 panels, may be in an ad- 
vanced post-ZAHB evolutionary state (see §7). The spectroscopic 
metallicities for these two stars currently are being firmed up and 
data reductions are continuing on the other 16 stars for which 
spectra are in hand. 
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Blazhko RR Lyr stars and the Lub-Simon galactic field stars (from 
S04) . The latter have metallicity estimates (on the ZW scale) from 
Layden (1994). The symbols are as indicated on the graph: Kepler 
stars (large black squares) ; two very metal-rich RR Lyr stars (red 
crosses); six metal-rich RR Lyr stars (red open squares); seven 
intermediate-metallicity RR Lyr stars (blue open triangles); and 
nine intermediate-to-low [Fe/H] RR Lyr stars (blue open circles). 
The Kepler stars have been shifted upwards by 0.151 (in accor- 
dance with Eq. 2) and thus are V values (as are the values for 
the Lub-Simon stars). 

most metal poor stars where the Fourier-based metallicity 
calibrations may give values systematically too high (see 
JK96 and Nemec 2004). The implication here is that the 
most metal-poor stars in our sample are even more metal 
deficient than the values given in Table 4. 

In general the metallicities that are based on the Fourier 
031 phase parameter (i.e., those derived using equation 1 of 
J98, and equation 3 of S04) are very similar. This can be 
seen by comparing the and (j!)§i values (ZW scale) given 
in columns 2 and 3 of Table 4. When the 031 metallicities 
are compared with those calculated using Atot (column 4) 
and RT (column 5) a number of systematic differences are 
present. For all the low-amplitude stars (the four suspected 
metal-rich stars, and AW Dra and KIC 7030715) the RT 
and j4tot estimates tend to be systematically more metal- 
poor than the 03i estimates. On the other hand, for most of 
the suspected metal-poor stars the Atot and RT equations 
give more metal-rich estimates of [Fe/H] than those derived 
with the 031 equations. And for several of the stars {e.g., 
V1510 Cyg, V346 Lyr, V350 Lyr, V894 Cyg and V1107 Cyg) 
the RT gives [Fe/H] values richer than the 03i values, while 
the Atot formula gives [Fe/H] values more metal poor. The 
largest difference occurs for V784 Cyg, where 03i suggests 
[Fe/H]zw = -0.16 while Aot and RT imply [Fe/H]zw = 
— 1 dex. Our suspicion is that the slope in the RT term of 
Equation 6 is too steep (see M3 line in lower-left panel of 
Fig.9). 

AW Dra is the only star for which the estimated [Fe/H] 
can be compared with a previously published metallicity es- 



timate, and, unfortunately, that estimate is quite uncertain. 
Based on its location in the P-Ay diagram, Castellani et al. 
(1998) found the metallicity of AW Dra to be in the range 
— 1.9 to —1.4 dex. Our estimate, [Fe/H]c9=— 1.67 dex, is 
consistent with theirs, but should be more reliable because 
it is based on much more extensive and higher precision 
photometry (see Table 3, and footnote 8). 

6.2 Reddenings and Mean Colours 

6.2.1 Dereddened colours 

J98 showed that the dereddened mean B — V colours corre- 
late well with pulsation period P and the Fourier Ai coeffi- 
cient, and derived the following relationship (her eq. 3): 

(B - 1/)o = 0.308 -f0.163P-0.187yli. (7) 

Replacing Ai{V) with Ai (Kp) from Eq. 2 simply reduces the 
intercept from 0.308 to 0.298. Dereddened B — V colours cal- 
culated with our modified equation for the ifp-photometry, 
and with the J98 equation for the l/-photometry of FN Lyr 
and AW Dra, are given in column 8 of Table 5 (under' J98'). 
The colour range is from 0.306 (for V368 Lyr) to 0.376 (for 
KIC 7030715), with a mean colour of 0.343. 

Another relationship for {B—V)q was derived by KWOl. 
Their eq. 6 includes an additive A3 term, and the Ai and A3 
are derived from ^-photometry. Using our offsets to convert 
to the Kp scale gives: 

[B - V)o = 0.448 + 0.189 logP - 0.313 Ai + 0.293 A3 

(8) 

where Ai and A3 now are derived from Kp photometry and 
the KWOl intercept has been modified accordingly. Dered- 
dened colours derived using this equation (and the KWOl 
equation for the V data of AW Dra and FN Lyr) are also 
given in column 8 of Table 5 (under 'KWOl'). These esti- 
mates are practically identical to the J98 colours and have 
an overall mean of 0.347. In subsequent analyses we adopt 
for each star the average of the two colour estimates. 

Figure 12 shows four graphs for the non-Blazhko 
RR Lyr stars, all with dereddened mean colour, [B — V)o, 
plotted along the x-axis. In every graph the four high 
metallicity stars (V782 Cyg, V784 Cyg, KIC 6100702 and 
V2470 Cyg), plotted as open squares, stand apart from the 
other RR Lyr stars. One sees that they also tend to be on the 
red side of the instability strip. The three stars with the red- 
dest colours (and coolest temperatures), NR Lyr, AW Dra 
and KIC 7030715, are all of very low metallicity and have the 
longest periods (P ~ 0.68 d) in our sample, and the three 
shortest period stars (V368 Lyr, V838 Cyg, V715 Cyg) have 
the bluest colours and are metal poor but of intermediate 
metallicity ([Fe/H]zw ~ —1.5 dex). Linear regressions were 
made to the points representing the 15 metal-poor stars, 
and the equations of the lines are given on the graphs. The 
colour-metallicity diagram (top left) shows that the reddest 
metal-poor stars tend to have the lowest metallicities; the 
colour-period diagram (top right) shows that the metal-poor 
stars exhibit a correlation such that the reddest stars have 
the longest periods, and that at a given colour the metal- 
rich stars have shorter periods than the metal-poor stars; the 
colour-total amplitude diagram (bottom left) shows that the 
red metal-poor stars have lower amplitudes than the blue 
metal-poor stars, and that at a given colour the metal-rich 
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Table 5. Galactic coordinates, reddenings, extinctions, dereddened colours and distances 



Star (Kp) {V) Z(II) 6(11) E(B - V) Ay (B - V)o d 

[mag] [mag] [dog] [deg] [mag] [mag] J98, KWOl [kpc] 



(1) 


(2) 


(3) 




(4) 




(5) 




(6) 


(7) 




(8) 




(9) 


NR Lyr 


12.683 


12.44 


69 


86686 


13 


55258 





175(0.009) 





532 





369, 





372 


3.0 


V715 Cyg 


16.265 




72 


92658 


07 


72890 





170(0.009) 





517 





321, 





327 


14.2 


V782 Cyg 


15.392 




75 


14435 


06 


94725 





509(0.024) 


1 


606 





358, 





363 


15.1 


V784 Cyg 


15.370 




76 


41126 


06 


56689 





343(0.010) 


1 


088 





351, 





353 


11.8 


KIC 6100702 


13.458 


13.64 


71 


0171 


17 


7007 





103(0.006) 





324 





348, 





353 


3.9 


NQ Lyr 


13.075 


13.36 


73 


11219 


15 


06421 





068(0.003) 





207 





352, 





358 


3.7 


FN Lyr 


12.876 


12.79 


73 


45862 


14 


69286 





104(0.007) 





320 





323, 





328 


2.9 


KIC 7021124 


13.550 




73 


5605 


14 


7208 





113(0.006) 





366 





357, 





361 


4.1 


KIC 7030715 


13.452 


13.24 


74 


5488 


12 


5345 





108(0.006) 





334 





376, 





377 


3.8 


V349 Lyr 


17.433 




72 


23472 


18 


37838 





102(0.013) 





301 





324, 





330 


22.7 


V368 Lyr 


16.002 




74 


40557 


14 


98618 





061(0.003) 





189 





306, 





308 


11.0 


V1510 Cyg 


14.494 




78 


7859 


07 


2303 





484(0.030) 


1 


480 





340, 





345 


10.4 


V346 Lyr 


16.421 




73 


65539 


19 


49070 





052(0.001) 





160 





340, 





345 


13.8 


V350 Lyr 


15.696 




75 


65435 


19 


58972 





051(0.001) 





158 





341, 





345 


9.9 


V894 Cyg 


13.393 


12.92 


78 


73702 


12 


60755 





112(0.005) 





346 





329, 





331 


3.2 


V2470 Cyg 


13.300 


13.54 


78 


3633 


14 


9140 





068(0.001) 





211 





357, 





360 


3.7 


V1107 Cyg 


15.648 




78 


54873 


15 


03295 





070(0.002) 





219 





348, 





355 


9.7 


V838 Cyg 


13.770 


14.24 


79 


20399 


16 


34690 





067(0.002) 





210 





314, 





318 


5.3 


AW Dra 


13.053 


12.85 


80 


23054 


19 


05392 





047(0.001) 





146 





363, 





365 


3.0 



stars have lower amplitudes than the metal poor stars; and 
the colour-absolute magnitude diagram (bottom right - see 
discussion in §6.4) shows that at a given colour the metal- 
rich stars have lower luminosity than the metal-poor stars, 
and that the reddest metal-poor stars are more luminous 
than the bluest metal-poor stars. 



6.2.2 Reddenings and extinctions 

Table 5 also contains E{B — V) reddenings (column 6) 
derived from the large-scale reddening maps of Schlegel, 
Finkbeiner & Davis (1998). These were estimated using 
the NStEd on-line facility and the galactic {l,h) coordinates 
given in Table 5 (columns 4-5). According to the maps the 
most reddened stars in our sample are V782 Cyg, V784 Cyg 
and V1510 Cyg. Total extinctions in the y-passband, Ay, 
were calculated assuming an extinction-to-reddening ratio 
Av/E{B — V) = 3.1 and are given in column 7 of Table 5. 

For AW Dra, which has calibrated B, V photometry, the 
above reddenings can be compared with reddenings derived 
from the Fourier-based dereddened colours and the observed 
reddened mean B — V colours. Castellani et al. (1998) ob- 
served (-B— y)o=0.368, and from the P-A diagram estimated 
[Fe/H] = -1.4 dex. Using the Burstein & Heiles (1982) red- 
dening maps they adopted the value E{B — V) = 0.04 — 0.06 
mag, which gives {B — l/)o=0.31-0.33. This estimate of the 
reddening is consistent with 0.047±0.001 given in Table 5, 
but the dereddened colour is bluer than the average of the 
two Fourier based values, {B - V)o= 0.364 (J98, KWOl). 
This difference explains why the Castellani et al. estimates 
of Tcff (see their table 7, and below) are higher than our 
estimates (column 2 of Table 6 below) . 



6.3 Effective Temperatures 

Mean T^s values for the RR Lyr stars were estimated using 
two prescriptions, one from KWOl and the other from S06 
(where the relationship between {B — V)o colour. Toff, and 
[Fe/H], including the effects of such other factors as surface 
gravity and turbulent velocity, has been reviewed). For the 
colour interval 0.20<(B - V)o<0.30 S06 adopts the formula 
given by Carney, Storm & Jones (1992): 

{B - V)o = -2.632 logToff + 0.038 [Fe/H]zw + 10.423. (9) 

Inverting this equation (eq. 18 of S06) gives: 

logToff = -0.380 {B~V)o + 0.0144 [Fe/H]zw + 3.960. (10) 

This equation, along with our colour and metallicity esti- 
mates, was used to calculate TeB for the 19 non-Blazhko 
stars. These are given in column 2 of Table 6 (under 'S06'). 
Since the Kepler RR Lyr stars have redder colours than 
the recommended range of applicability of the S06 equation 
Tcff values also were calculated using eq. 11 of KWOl. For 
[«/Fe]=0 the metallicity term [M/H]cg is equal to [Fe/H]cG 
and we can write 

logTeff = 3.8840 - 0.3219 {B - V)o + 0.0167 log 5 
-F0.0070 [Fe/H]cG, 

which is valid over the colour range of the ab-type RR Lyr 
stars studied here. Estimates of Tcff based on this formula 
also are given in Table 6 (column 2, under KWOl). In general 
the temperatures agree, with an average difference of only 
20 K, and one could simply take the average for the best es- 
timate; however, the KWOl temperatures are higher (lower) 
than the S06 temperatures for the metal-poor (metal-rich) 
stars. We shall see below that there is reason to believe that 
the KWOl formula is to be preferred (see §6.9, and in par- 
ticular the bottom left panel of Fig. 13). 

Castellani et al. (1998, Table 7) found Tcff ~ 6700 K 
for AW Dra, assuming fundamental mode pulsation, mass 
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Figure 12. Four diagrams for the non-Blazhko RR Lyr stars, all with dereddened mean colour {B — V)o (average of the two values 
given in column 8 of Table 5) along the abscissa: metallicity vs. colour (Top left), period vs. colour (Top right), total amplitude vs. colour 
(Bottom left), and absolute magnitude vs. colour (Bottom right). In each panel the line (and its equation) is from a least squares fit 
to the points for the 15 metal-poor stars (i.e., those with [Fe/H]< —1.0 dex). The linear correlations show that the reddest metal-poor 
stars (those nearest the red edge of the Instability Strip) are more metal poor, have longer periods and smaller total amplitudes, and are 
more luminous than the bluest metal-poor stars. The diagrams also show that at a given colour the four metal-rich stars (plotted with 
open squares) have shorter periods, smaller amplitudes and are less luminous than the metal-poor stars. 



~0. 65-0. 75 Mq, and luminosity ~65-80 Lq. This tempera- 
ture is ~400 K hotter than the corresponding Fourier-based 
estimate of Tea ~ 6300 K. The difference is not surprising 
given that their assumed mass and luminosity are higher 
than the Fourier-based values. 



6.4 Absolute magnitudes 

Following on from the Jurcsik & Kovacs (1996) and Kovacs 
& Jurcsik (1996, 1997) papers, J98 gives an equation for 
absolute magnitude (based on V-photometry) that depends 
on Ai and (her eq. 2). Application of our V-Kp offsets 
(Eq. 2) gives 

Mv = 1.179 - 1.396 P- 0.477 Ai +0.103 </>ii, (12) 



where Ai and (fif^i are Fourier parameters derived from the 
.R'p-photometry. The resulting Mv values are given in col- 
umn 3 of Table 6 (under 'J98, eq2'). 

The Fernley et al. (1998) equation relating Mv and 
metallicity, 

Mv = (0.20 ± 0.04) [Fe/H]cG + (1-03 ± 0.14), (13) 

assumes that RR Lyr has Mv = 0.78(±0.29) at [Fe/H]cG = 
— 1.39 dex and is consistent with the statistical parallax so- 
lution for 84 halo RR Lyr stars and with various Baade- 
Wesselink analyses. Mv values derived using this equation 
are listed in column 3 of Table 6 (under 'F98, eq.2'). In 
general the F98 and J98 values are similar, with a mean 
difference (J98 minus F98) of only 0.01(±0.01) mag. On the 
downside there has been mounting evidence that the calibra- 
tion of Mv with [Fe/H] is non-hnear (see S06 and TSR8). 
Equation 10 of Bono, Caputo & Di Criscienzo (2007) 
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Table 6. Fundamental physical properties of the Kepler non-Blazhko RR Lyr stars 



Star 




Mv 


log 9 


L(puls) 


L(evol) 


A^(puls) 


A^(evol) 




KWOl, S06 


J98, F98, BCD7, CCS 


J98, KW99 


J98, J98 


S06, S06, S06 


J98, J98 


S06, BCD7 




eqll, eql8 


eq2, eq2, eqlO, eq4a 


eql5, eql2 


eql6, eql7 


eq8, eqlO, eql2 


eql4, eq22 


eql5, eq7 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



NR Lyr 


6242 


6124 





63 





64 





67 





47 


2 


677 


2 


689 


46 


9 


47 


3 


58 


7 


54 


I 


61 


I 





66 





62 





73 





72 


V71 ^ Cw 


6607 


6603 


Q 


86 


Q 


81 


Q 


81 


Q 


67 


2 


874 


2 


882 


38 


Q 


36 




47 


g 


45 


\ 


48 


Q 


Q 


60 


Q 


55 


Q 


64 


Q 


64 




6493 


6512 





96 





95 


I 


03 





84 


2 


817 


2 


826 


32 





34 


g 


30 


5 


34 


4 


40 


3 





53 





50 





57 





58 


V784 Cyg 


6546, 


6577 





95, 





98, 


1 


08, 





88 


2 


807, 


2 


815 


31 


0, 


33 


3 


27 


2, 


32 


2, 


38 


9 





49, 





49 





56, 





56 


KIC 6100702 


6564, 


6605 





99, 





98, 


1 


08, 





87 


2 


855, 


2 


863 


31 


0, 


33 


1 


27 


4, 


32 


3, 


38 


9 





52, 





49 





56, 





56 


NQ Lyr 


6382, 


6358 





75, 





73, 





72, 





57 


2 


756, 


2 


767 


42 


1, 


42 


2 


55 


2, 


50 


3, 


54 


3 





63, 





58 





68, 





68 


FN Lyr (Kp) 


6531, 


6482 





76, 





72, 





71, 





56 


2 


814, 


2 


822 


42 


7, 


40 


2 


55 


8, 


50 


8, 


55 


1 





60, 





59 





69, 





68 


FN Lyr (V) 


6559, 


6551 





76, 





72, 





72, 





56 


2 


814, 


2 


823 


42 


5, 


39 


4 


55 


6, 


50 


6, 


54 


8 





60, 





59 





69, 





68 


KIC 7021124 


6332, 


6240 





70, 





69, 





69, 





52 


2 


725, 


2 


737 


44 


5, 


44 


4 


57 


6, 


52 


4, 


57 


6 





64, 





60 





71, 





70 


KIC 7030715 


6267, 


6181 





69, 





74, 





73, 





58 


2 


676, 


2 


688 


41 


7, 


43 


7 


54 


5, 


49 


8, 


53 


7 





58, 





58 





68, 





67 


V349 Lyr 


6550, 


6517 





80, 





76, 





74, 





60 


2 


835, 


2 


843 


40 


8, 


39 





53 


2, 


48 


8, 


52 


5 





61, 





57 





67, 





66 


V368 Lyr 


6683, 


6689 





84, 





80, 





79, 





65 


2 


891, 


2 


898 


38 


9, 


36 


2 


49 


4, 


46 


3, 


49 


8 





57, 





56 





65, 





65 


V1510 Cyg 


6435, 


6370 





73, 





72, 





71, 





56 


2 


762, 


2 


772 


42 


8, 


41 


8 


56 


0, 


51 


0, 


55 


3 





60, 





59 





69, 





68 


V346 Lyr 


6439, 


6376 





74, 





72, 





71, 





56 


2 


766, 


2 


776 


42 


6, 


41 


6 


55 


8, 


50 


8, 


55 








60, 





59 





69, 





68 


V350 Lyr 


6426, 


6359 





72, 





71, 





71, 





55 


2 


750, 


2 


760 


43 


1, 


42 


1 


56 


3, 


51 


2, 


55 


6 





58, 





59 





69, 





69 


V894 Cyg 


6496, 


6443 





72, 





71, 





71, 





55 


2 


771, 


2 


781 


42 


9, 


41 





56 


2, 


51 


1, 


55 


5 





56, 





59 





69, 





68 


V2470 Cyg 


6469, 


6470 





90, 





91, 





94, 





79 


2 


793, 


2 


802 


33 


9, 


36 


1 


36 


4, 


38 


1, 


42 


9 





53, 





51 





59, 





60 


V1107 Cyg 


6425, 


6375 





79, 





77, 





76, 





62 


2 


776, 


2 


786 


40 


1, 


40 


3 


51 


8, 


47 


9, 


51 


4 





61, 





57 





66, 





66 


V838 Cyg 


6627, 


6617 





82, 





78, 





77, 





64 


2 


863, 


2 


871 


39 


5, 


37 


3 


50 


8, 


47 


2, 


50 


7 





58, 





56 





66, 





65 


AW Dra (Kp) 


6306, 


6217 





65, 





71, 





70, 





54 


2 


673, 


2 


685 


43 


4, 


44 


1 


56 


7, 


51 


5, 


56 


1 





55, 





59 





70, 





69 


AW Dra (V) 


6303, 


6276 





65, 





70, 





70, 





54 


2 


673, 


2 


684 


43 


3, 


43 


3 


56 


6, 


51 


4, 


56 








51, 





59 





70, 





69 



provides a third estimate of Mv- They show that over the 
entire metallicity range from [Fe/H]zw = —2.5 to 0.0 dex 
the galactic field variables with Layden (2007) metallicities 
are described by 

Mv = 1.19 (±0.10) + 0.50 [Fe/H]2w + 0.09 [Fe/HJ^w (14) 

Absolute magnitudes calculated with this quadratic formula 
also are given in column 3 of Table 6 (under 'BCD7'). For the 
metal-poor stars the BCD7 A4v values are very similar to the 
F98 values; however, because of the non-linear metallicity 
dependence they are fainter than the F98 values for the four 
metal-rich stars. 

Most recently, Catelan & Cortes (2008) argue that re- 
vised values for the trigonometric parallax and reddening of 
RR Lyr imply that the luminosity scale for RR Lyr stars 
should be brighter and give the following equation: 

Mv = (0.23 ± 0.04) [Fe/H]zw + (0.984 ± 0.127). (15) 

Absolute magnitudes calculated with this equation (given in 
column 3 of Table 6 under 'CC8') are, on average, ~0.15 mag 
brighter than the values from F98 and BCD7, which in turn 
are ~0.05 brighter than the J98 values. Note that these dif- 
ferences are smaller than the uncertainties in the individual 
Mv estimates (which are ~0.20 mag). The colour- Mv dia- 
gram plotted in Fig. 12 (bottom right panel) was constructed 
using the Mv values calculated with eq. 2 of J98. 

6.5 Distances 

Approximate distances for the non-Blazhko RR Lyr stars 
are given in column 9 of Table 5. These were computed as- 
suming: (1) the Fernley et al. (1998) Afv values (which lie 



between the J98 and CC08 values); (2) the AS AS (V) val- 
ues (given in column 3 of Table 5) for the nine stars having 
ASAS V, I photometry; and for the other stars the (Kp)- 
magnitudes brightened by 0.15 mag (the observed average 
V-Kp offset for FN Lyr and AW Dra); and (3) the visual 
extinctions given in column 7 of Table 5. With uncertain- 
ties of ~0.15 mag for Mv, ~0.08 mag for the V-Kp offset, 
and ~0.03 mag for Av, these distances are quite uncertain. 
Taken at face value the nearest of the stars have distances 
~3 kpc, and the most distant star is V349 Lyr at ~22 kpc. 
These estimates will be considerably improved when cali- 
brated _B\/ii/- photometry becomes available. 

6.6 Surface gravities 

Mean surface gravities were calculated using eq. 15 from J98, 
and eq. 12 from KW99. The J98 formula, which is accurate 
to ±0.004, depends only on the period and is given by 

log g = 2.473 - 1.226 logP. (16) 

The KW99 formula depends on period, mass and effective 
temperature: 

log.g = 2.938 + 0.230 logTM - 0.110 logTcfi - 1.219 logP 

(17) 

Both estimates are given in column 4 of Table 6 and the 
results are in good agreement. 

6.7 Pulsational luminosities and masses 

Following on from the basic equation of stellar pulsation 
given by the Ritter (1879) relation, Py^ — Q (where p is 



Kepler non-Blazhko RR Lyr Stars 27 




3.90 

0.06 
0.04 
0.02 
0.00 
-0.02 
-0.04 
-0.06 
-0.08 
-0.10 
-0.12 



M: J98, eq. (22); L: J98 eq. (17) 



3.80 



0.45 



0.55 0.6 
P [days] 



0.65 



3.75 










• 




• 






• 














o • 

•• 




• 






• 









So 


• 


• 


• 


• 


• - 

o - 


9 • 


G 




o o 








-KWOl • 
S06 O 








O 












o ■ 



0.7 




[Fe/Hl-0.47, [O/Fe|=0.23, Y=0.266 
[Fe/H]=-0.47, [O/Fel=0.0, Y=0.266 
[Fe/H]-1.7B, [O/Fe]=0.66, Y=0.24B 
[Fe/H]-1.7B, [O/Fel=0.0. Y-0.247 



3.90 



0.02 
0.00 
-0.02 
-0.04 
-0.06 
-0.08 
-0.10 
-0.12 
-0.14 
-0.16 
-0.18 
-0.20 



3-85 log T,„ 

M: S06, eq. (15); L: S06 eq. (12) 



3.80 



3.75 



1^ 


• 






• - 








• 




2 • 






O - 


• 











• • 








• • 




• 




• 


• 












• 

8 • 








% m. ,n o 




o 


. o 



0.45 



0.5 



0.55 0.6 
P [days] 



Figure 13. (Top panels) HR-diagrams showing tlic locations of tlic Kepler non-Blazhko RR Lyr stars compared with two sets of 
horizontal branch models (low- and high-metallicity) from Dorman (1992). In both graphs the four metal-rich RR Lyr stars have been 
plotted with open squares and the metal-poor stars with filled squares; the black numbers next to the ZAHB symbols and the two 
evolutionary tracks are the assumed masses for the individual stellar evolution models. For the non-Blazhko RR Lyr stars in the top left 
panel the L(puls) were calculated with eq. 17 of J98 and labelled (in red) with the A^(puls) calculated with eq. 22 of J98; in the top 
right panel the L(evol) were calculated with eq. 12 of S06 and labelled (in red) with the jV((evol) calculated with eq. 15 of S06. Both 
panels (left and right) show the blue and red edges for the instability strip computed using the Warsaw pulsation code; the fundamental 
mode edges are plotted as solid lines, and first-overtone mode edges as dashed lines. (Bottom panels) Graphs comparing the observed 
pulsation periods and the pulsation periods derived with the Warsaw pulsation code. The assumed masses and luminosities are the values 
given in the top panels. For each star two points are plotted, computed assuming the Tgfj from eq. 11 of KWOl (blue dots) and from 
eq. 18 of S06 (red circles). The best agreement is seen in the left panel, with the KWOl temperatures favoured over the S06 values. In 
both panels the largest differences are seen for the longest period stars, with no metallicity dependence. 



the mean density in cgs units and Q is the pulsation con- 
stant), van Albada & Baker (1971, 1973) derived an equation 
relating the pulsation period to mass, luminosity and T^s. 
When this equation and similar more recent equations that 
include a dependence on metal abundance are used to derive 
the mass and luminosity for an RR Lyr star such quantities 
are referred to as pulsational mass M (puis) and pulsational 
luminosity L(puls). Both are expressed here in solar units. 

Pulsational luminosities were calculated with two differ- 
ent formulae given by J98 and are reported in column 5 of 
Table 6. In the first case (her eq. 16) the luminosity depends 
only on the metallicity; 

logL = 1.464 - 0.106 [Fe/H]cG. (18) 

The second formula (her eq. 17) also takes into account T^a'- 

logL = 10.260 - 0.062 [Fe/H]cG - 2.294 logTes. (19) 

In both cases the lower metallicity stars have the higher lu- 
minosities. In the second equation the KWOl effective tem- 



peratures (Table 6, column 2) were used for the calculations. 
With both equations one sees in the HR-diagram two ap- 
proximately parallel lines, one for the metal-rich stars and 
one for the metal-poor stars, each with the luminosity in- 
creasing as Toff decreases (see the top left panel of Fig. 13). 
At a given T^s the luminosity difference between the metal- 
rich and metal-poor stars is greater with eq. 16 than it is 
with eq. 17. 

Column 7 of Table 6 contains M (puis) values computed 
using eq. 14 and 22 of J98. The former is given by 

logX = 1.477 logL - 1.754 logP - 6.272 logT.ff 
+ 0.037 [Fe/H]cG + 20.884, 

and the latter by 

logX = -0.328 - 0.062 [Fe/H]cG. (21) 

In both cases the adopted [Fe/H] values given in Table 4 
(transformed to the CG system) were used, and in the first 
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equation we used the logL from J98 eq. 17 and the logTcff 
from eq. 11 of KWOl. The average mass for the four metal- 
rich stars is ~0.50 Mq compared with the average mass for 
the metal-poor stars of ~0.60 Mq- 



6.8 L and A4 from stellar evolution models 

A ZAHB stellar evolution model takes as its input the mass, 
A^(evol), and chemical composition {X, Y,Z). A subsequent 
evolutionary track for a given mass and composition gives 
the luminosity L(evol) and effective temperature as a func- 
tion of time. Examples of such models are those by Dorman 
(1992), Bono et al. (1997) and VandenBerg et al. (2000). 
Based on three different sets of stellar evolution models S6 
gives formulae for L(evol) as a function of [Fe/H]. His eq. 8, 
which follows from the models of Caputo et al. (2000), is 
given by 

logL = 1.245 - 0.451 [Fe/H] - 0.097 [Fe/Uf. (22) 

His eq. 10, which is based on the alpha-enhanced ZAHB 
models of Catelan, Pritzl & Smith (2004), is given by 

logL = 1.404 - 0.243 [Fe/H] - 0.043 [Fe/H]^ (23) 

Both of these equations, like the BCD7 absolute magnitude 
formula discussed above, assume a quadratic dependence on 
[Fe/H]. This is not the case for eq. 12 which is derived from 
the models of Clementini et al. (2003) and is given by 

logL = 1.538 - 1.110 [Fe/H]. (24) 

Luminosities computed with these three formulae are given 
in column 6 of Table 6. For the metal-poor stars all three 
L(evol) are systematically larger than the L(puls) values 
given in column 5. For the four metal-rich stars the agree- 
ment is better but there is a wide range of L(evol) owing to 
the linear or non-linear [Fe/H] dependencies. Regardless of 
which formula was used the most luminous stars have the 
lowest metallicities and there is internal consistency (as was 
the case for the J98 luminosities). It is not clear whether 
the L(puls) or L(evol) are correct. Further comparison of 
the derived luminosities is given in §6.9 (after masses are 
discussed). 

A mass equation with A4 varying linearly with [Fc/H] 
is also given by S06 (his eq. 15): 

logX = -0.283 - 0.066 [Fe/H]zw. (25) 

Since this equation was derived from the Bono et al. (1997) 
horizontal branch models it gives evolutionary masses. The 
Al(evol) derived with this equation are given in column 8 
of Table 6. A more recent mass formula by Bono, Caputo 
& Di Criscienzo (2007), their eq. 7, which is based on the 
Pietrinferni et al. (2004, 2006) HB models, is given by 

(logM) = -0.2675 - 0.063 [Fe/H]. (26) 

This formula is very similar to the S06 equation. The masses 
computed with this formula are given in column 8 of Ta- 
ble 6 and are seen to agree to within 0.01 Mq with the S06 
masses. As was the case for the luminosities, the M{evol) 
are all larger than the corresponding Al(puls) values. 



6.9 Comparison with evolutionary models 

Having computed M , L and T^b values for the non-Blazhko 
RR Lyr stars from observed Fourier parameters using equa- 
tions that derive from both pulsation and stellar evolution 
theory it is of interest to compare the results with the loca- 
tions of model ZAHBs and evolutionary tracks. 

The top two panels of Figure 13 show HR diagrams 
with logTefi as abscissa and logL as ordinate, and two sets 
of ZAHB loci from Dorman (1992) computed for differ- 
ent masses along the ZAHB. The more luminous horizontal 
branch assumes [Fe/H]= —1.78 dex, and the less luminous 
branch [Fe/H]= —0.47 dex. The numbers next to the sym- 
bols are the masses for the individual models, which are 
seen to be higher for the low-[Fe/H] tracks than for the 
high-[Fe/H] tracks. For both assumed metallicities oxygen 
enhanced and non-enhanced ZAHBs have been plotted - the 
effect of increasing the oxygen to iron ratio from [O /Fe]=0 to 
0.66 for the low-metallicity ZAHBs is to lower the luminosity 
and reduce the mass at a given temperature. For the high- 
metallicity tracks an oxygen enhancement from [O/Fe]=0 to 
0.23 has little effect on the derived L or M. Also plotted in 
the low- [Fe/H] case are the evolutionary paths away from 
the ZAHB for two masses, 0.66 and 0.68 Mq. In both pan- 
els the non-Blazhko RR Lyr stars with low metallicities are 
represented by large black squares, the four high-metallicity 
stars are plotted with open squares, and the T^s are the 
average of the KWOl and S06 values (column 2 of Table 6). 

In the top left panel of Fig. 13 the luminosities and 
masses (labelled in red) of the Kepler RR Lyr stars were 
calculated with eqs. 17 and 22 of J98 and thus are based on 
pulsation theory. The L(puls) and Al(puls) are seen to be 
systematically smaller than values derived from the ZAHB 
tracks (for the appropriate metal abundance). The reddest 
non-Blazhko RR Lyr stars lie close to the fundamental mode 
red-edge and have the smallest amplitudes. This graph also 
shows blue and red edges of the instability strip for the fun- 
damental mode (red and blue solid lines) and first-overtone 
mode (red and blue dashed lines). The edges were calculated 
with the Warsaw pulsation code (see Section 7) assuming a 
mass of 0.65 Mq. The Kepler non-Blazhko stars all lie in 
the fundamental mode region of the variability strip, and the 
smallest amplitude RR Lyr stars (the four metal-rich stars, 
KIC 7030715 and NR Lyr) have locations near the funda- 
mental red edge (FRE) of the instability strip. As expected, 
all the stars near the FRE have low Rsi values. 

In the top right panel of Fig. 13 the luminosities were 
calculated with eq. 12 of S06, and the masses (labelled in 
red) with eq. 15 of S06; thus they are evolutionary L and 
M values. In this case there is very good agreement with 
the stellar evolution models, as one expects since they are 
based on stellar evolution models. Enhancing the oxygen to 
iron ratio by the plotted amounts makes little difference. 

It is unclear which are correct, the L(puls) and 
A/f(puls), or the L(evol) and Al(evol)? The mass and lu- 
minosity discrepencies go in the same direction as seen for 
Cepheids. Pietrzynski et al. (2010) recently derived a dy- 
namical mass 7Vf(dynam) for a classical Cepheid in a well 
detached, double- lined eclipsing binary in the LMC. The 
mass they derive is very accurate and favours M (puis) . The 
reason for the discrepancies may be the same, as suggested 
by Pietrzynski et al. - not enough mass loss has been taken 
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into account in the evolution models. For further guidance 
on these questions we turn to pulsation models. 



7 HYDRODYNAMIC MODELS 

Smolec & Moskalik (2008) recently have developed the War- 
saw convective pulsation programs for studying stellar pul- 
sation. The codes, both linear and non-linear, are one di- 
mensional and use a single equation to describe the genera- 
tion of turbulent energy; this is done according to the model 
proposed by Kuhfufi (1986, see also Wuchterl & Feuchtinger 
1998). Even though a simple diffusion approximation is used 
to describe the radiation field the models are able to repro- 
duce quite well the dynamics of RR Lyrae pulsations. 

The most recent application of these programs (Smolec 
et al. 2011) has been to construct hydrodynamic models for 
the purpose of testing Stothers' (2006) proposed explana- 
tion of the Blazhko phenomenon that is observed in about 
half of all RR Lyr stars. In these models the strength of the 
turbulent convection was modulated and the resulting mod- 
els were compared in detail with the Fourier descriptions of 
the Kepler observations of the RR Lyr stars (Kolenberg et 
al. 2011). 

The initial unmodulated models of the Smolec et al. 
(2011) study are also applicable for the non-Blazhko stars 
considered here. The predicted bolometric light curves were 
transformed to l/-band light curves by applying to each pul- 
sation phase a bolometric correction derived from Kurucz 
(2005). 

7.1 Pulsational or Evolutionary L, M? 

As a check on the luminosities and masses derived from stel- 
lar pulsation theory and from the horizontal branch mod- 
els the Warsaw pulsation code was used to calculate pulsa- 
tion periods. Such periods are very robust and depend very 
weakly on the assumed convective parameters. The calcula- 
tions assumed the metallicities given in Table 4, the L(puls) 
and Al(puls) values (calculated with the J98 equations) 
given in the top left panel of Fig. 13, and the L(evol) and 
M{evol) values (from the S06 equations) given in the top 
right panel of Fig. 13. The results, summarized in 'percent- 
difference' diagrams comparing the observed pulsation peri- 
ods and the model pulsation periods, are shown in the bot- 
tom panels of Fig. 13. Because the assumed effective temper- 
atures are important, in both panels two points are plotted 
for each star, one computed assuming the Tcb from eq. 11 
of KWOl (blue dots), and the other assuming the Teg from 
eq. 18 of S06 (red circles). For the pulsational values (left 
panel) the best formulae seem to be those in J98, and the 
KWOl temperatures (blue dots) seem to be preferred over 
the S06 temperatures (red circles) - in this case the observed 
and calculated periods always agree to within 5% and ~half 
the stars are within 2%. On the other hand, the L and A4 
values from the stellar evolution models (right panel) result 
in model periods systematically too large by as much as 20%, 
with almost no agreement. In both panels the largest differ- 
ences between the KWOl and S06 temperatures are seen 
for the longest period stars. The four metal-rich stars all 
have periods shorter than 0.55 d and there is no systematic 
metallicity dependence. 



7.2 Predicted Fourier parameters 

Three sets of models were computed with the Warsaw code, 
varying L while holding the other variables constant, then 
M, then [Fe/H]. The resulting Fourier parameters were 
plotted in diagrams such as those found in Figs. 2 and 
6-8. In Figure 14 the luminosity was varied {L/Lq = 
40, 50, 60) while the metallicity was held constant at [Fe/H] = 
— 1.30 dex {Z — 0.001) with the hydrogen fraction fixed 
at X=0.76. Two masses were considered, M/Mq = 0.65 
(solid lines) and 0.55 (dashed lines). In Figure 15 the 
metallicity was varied from [Fe/H]= —0.24 dex {Z = 0.01), 
to -1.24 dex (Z = 0.001), to -1.93 dex {Z = 0.0002)0 
while keeping the other variables constant at L/Lq — 50, 
M/Mq = 0.65, and X=0.76. 

The Al vs. yltot diagrams (top left panels of Figs. 14-15) 
show very good agreement with the estimates for the non- 
Blazhko stars (right panel of Fig. 2). The approximately 
linear trend seems to be independent of luminosity or mass 
or metallicity variations. However, the fit is not exact and 
the observed non-linearities (which are small) appear to go 
in opposite directions to the model predictions. 

The 4)21 vs Al graphs (middle left panels of Figs. 14- 
15) appear to be quite sensitive to L, M and metallicity 
effects (see the lower left panels of Figs. 6-7). A star of given 
mass and composition will move up and to the right as its 
luminosity increases. Alternatively, a star of given L and 
composition will move up and to the right if the mass de- 
creases. And, if L, M and X are kept constant metal-rich 
stars ought to occupy the high Ai, high (j!)2i region of the 
diagram. 

To attempt an application consider the stars AW Dra 
and V784 Cyg, both of which are seen in Fig. 6 with higher 
than average 4)2i values at a given Ai. If the models are 
indicative, then the two stars with the largest luminosities 
might be expected to be AW Dra and V784 Cyg. According 
to column 5 of Table 6 AW Dra is one of the most luminous 
stars in the sample, but V784 Cyg is not - it has one of 
the lowest luminosities. On the other hand, V784 Cyg is one 
of the four metal-rich stars, suspected of having low L and 
low mass. The high location in Fig. 6 of V784 Cyg (and the 
other high metallicity stars) might be explained by its low 
mass, whereas the location of AW Dra would seem to be due 
primarily to its high luminosity (either as a low- [Fe/H] star 
near the ZAHB or as a more metal-rich star in an advanced 
evolutionary state). The highest L (and highest M) star 
in the sample, NR Lyr (see Table 6), also has the lowest 
metallicity (see Table 4), even lower than that of AW Dra; 
for it the lower [Fe/H] would seem to cancel out the higher L, 
thus explaining its smaller amplitude and smaller 021 than 
AW Dra. 

The top right and middle right diagrams of Figs. 14- 15 
are somewhat noisy. The top-right panel of Fig. 14 suggests 
that the lowest L stars have the smallest R21 values (note 
that the scale matches that shown in the upper right panel 
of Fig. 6, but is much reduced from that shown in Fig. 7 for 
globular cluster stars). Basically, R21 and Rsi have to be 
small very close to the red edge of the instability strip. At 

^ When transforming between Z and [Fe/H] we used [Fe/H] = 
logZ — logZQ = logZ + 1.765, where the solar metallicity was 
assumed to be Zq = 0.01716 (Sweigart & Catelan 1998). 
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the red edge the amphtude goes to zero, as do R21 and -R31. 
The middle right panels of Fig. 14 suggest that the lowest-L 
stars are not expected to have high 02i values. This is to 
be compared with the Kepler data shown in the lower right 
panel of Fig. 6 that suggest that there is a nearly linear 
relationship between (p2i and 03i - this is not seen in the 
models. 

The bottom left panels of Figs. 14-15 show that at a 
given 031 a shift to the right can be caused by a luminosity 
increase or a higher mass. The bottom left panel of Fig. 15 
shows the traditional result that at a given 031 metallicity 
decreases as period increases; but it also seems to be sug- 



gesting that at a given period 03i increases with decreasing 
metallicity. 

Finally, the bottom right panel in Fig. 14 shows that 
at a given yltot and metal abundance stars of higher L for 
a given mass, or of lower mass for a given L, have longer 
periods; this is as expected, and similar results are seen in fig. 
6 of DairOra et al. (2003). The bottom right panel of Fig. 15 
shows that for a given L and A4 varying the metallicity has 
a relatively small effect on the period. Taken together these 
two panels suggest that the main factor shifting the periods 
to longer values in the Bailey diagram (or any of its surrogate 
diagrams - see Fig. 9) seems to be higher luminosities (which 
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Figure 15. Graphs showing the effect on the Fourier parameters derived with the Warsaw convective pulsation code of varying [Fe/H] 
from -0.24 dex (Z = 0.01, green), to -1.24 dex (Z = 0.001, blue), to -1.93 dex (Z = 0.0002, red), while keeping the mass, luminosity 
and hydrogen content constant at 0.65Mq, 50I/q and X=0.76. All ordinates are for the V-passband. 



can be caused either by higher ZAHB mass or by advanced 
HB evolution) and not lower metallicities. 

Obviously more work needs to be done to optimize the 
estimation of the physical variables in these observational 
planes, but the potential seems high. 



8 SUMMARY 

The main results of this paper are as follows: 

(1) Fourier decomposition has been performed on the 
19 least modulated ab-type RR Lyr stars observed with the 
Kepler space telescope at long cadence (every 30 min) and 



short cadence (every 1 minute) during the first 417 days of 
its operation (Q0-Q5); 

(2) While none of the RRab stars shows the recently 
discovered 'period-doubling' effect seen in Blazhko variables, 
the star KIC 7021124 was discovered to pulsate in the fun- 
damental and second overtone modes with a period ratio 
P2/P0 ~ 0.59305 and to have properties similar to those of 
V350 Lyr; 

(3) Period change rates and improved periods have been 
derived from O— C diagrams for several of the stars that have 
historical data. For AW Dra, data from the last 12 years 
suggest that its period is increasing at the rate dP/d£=3.79 
d/Myr, while data spanning the last 100 years suggest a 



32 Nemec, Smolec, Benko et al. 



slower rate, dP/dt=0.32 d/Myr. FN Lyr appears to have a 
very slowly increasing period, and the periods of NR Lyr 
and NQ Lyr appear to be constant. 

(4) Because the differences between the Kp and V 
Fourier parameters for three stars (AW Dra, FN Lyr an 
NR Lyr) are found to be small and systematically different 
we were able to use extant l/-band correlations (with small 
modifications) to derive underlying physical characteristics 
for the Kepler stars. This procedure seems to be validated 
through comparisons of the Kepler stars with other galac- 
tic and LMC field RR Lyr stars and with RR Lyr stars in 
galactic and LMC globular clusters. 

(5) Preliminary metal abundances have been derived for 
all the non-Blazhko RR Lyr stars. Thirteen of the stars ap- 
pear to be similar to those found in intermediate metallicity 
globular clusters (i.e., [Fe/H]c9~ —1.6 dex); the most metal 
deficient star appears to be NR Lyr with [Fe/H]c9 ~ —2.3 
dex, and the four lowest amplitude stars (KIC 6100702, 
KIC 9947026, V782 Cyg and V784 Cyg) appear to be metal- 
rich with [Fe/H]c9 between —0.55 and +0.07 dex. 

(6) In general the luminosities of the metal-rich RR Lyr 
stars are found to be lower than those of the metal-poor 
stars; however, the luminosities derived from stellar evolu- 
tion models are systematically higher than those derived 
from stellar pulsation models. It is not clear which are cor- 
rect. 

(7) The three stars with the longest periods, AW Dra, 
NR Lyr and KIC 7030715 (ah with periods ~0.68 d) also 
are the reddest and coolest stars. We suspect that AW Dra 
may be in an evolved state. In general, the stars with the 
lowest amplitudes are found to be located nearest the red 
edge of the instability strip. 

(8) The mass range for the entire sample (approximate) 
is from 0.50 to 0.65 Mq if based on pulsation theory, or from 
0.56 to 0.72 if based on HE evolution models. 

(9) Finally, the Fourier parameters of the stars have 
been compared with values newly computed with the War- 
saw convective pulsation codes. We find that in the P-A 
diagram varying the metallicity for a given L and A4 has a 
relatively small effect on the period-shift at a given ampli- 
tude (in fact, the small shift is toward shorter periods for 
more metal poor stars), and that the main factors causing 
the period shifts must be L and A4. 
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SUPPORTING INFORMATION 

To illustrate the constancy of the pulsations of the stars 
studied here 'animated gifs' have been prepared for all the 
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stars. These may be found in the onhne version of this arti- 
cle, along with the data files. 



APPENDIX A: KIC 7021124 

During the course of our analyses the star KIC 7021124 
was discovered to be a doubly-periodic RR Lyr star, sim- 
ilar to V350 Lyr (see BIO). From an analysis of the Ql 
data (1626 long cadence observations made over a 30-day 
period) its period was found to be Pq = 0.6224925(7) d, 
corresponding to the frequency fo = 1.606445(2) d^^. A 
light curve phased with this period is plotted in the upper 
panel of Figure Al. The bottom panel shows the Fourier 
transform after prewhitening with fo and its harmonics. In 
addition to residual power seen at the location of the re- 
moved primary frequency and its harmonics, we also see a 
family of at least six additional peaks (labelled with arrows). 
Adopting the feature at 2.7088 d~^ as the independent fre- 
quency /2, the other frequencies correspond to /2 ± fc/o, 
where k is an integer. With this assumption the period ratio 
is P2/P0 = 0.59305, which is almost identical to that derived 
by BIO for V350 Lyr {P2/P0 = 0.592). The two stars also 
are similar in their fundamental periods, in their Fourier 
characteristics (see Figs. 6,9,10), and in their masses and 
luminosities (see Table 5). 

Figure A2 shows KIC 7021124 (and V350 Lyr) plot- 
ted in the P2/P0 vs. Po 'Petersen' diagram. The curves 
were computed with the Warsaw pulsation hydrocode in- 
cluding turbulent convection. The periods and period ratios 
do not depend strongly on the convective parameters en- 
tering the model (which in this case were set C, adopted 
in Baranowski et al. 2009). Other model parameters are: 
Z=0.0001 (or [Fe/H]=-2.2); X=0.76 (latest opacities, with 
new solar mixture); three masses - O.55A10 (open circles), 
0.65Mq (filled circles), 0.75Mq (crosses); and four lumi- 
nosities: 401/0 (red), 501/0 (green), 60Lq (blue), 701/0 (pur- 
ple). The best agreement for KIC 7021124 (and V350 Lyr) is 
obtained for a high luminosity and a high mass: L/Lq = 70 
and M/Mq = 0.75. In Table 6 both stars are among the 
highest L and M stars in our sample. The high values in- 
ferred from the Petersen diagram are more in accord with 
the evolutionary values than the pulsation values. 



APPENDIX B: ASAS l^, /-PHOTOMETRY 

Very little colour information is available for the non- 
Blazhko RR Lyr stars. Fortunately the ASAS-North sur- 
vey (see Pigulski et al. 2009) includes calibrated V and / 
photometry for nine of the brightest Kepler non-Blazhko 
stars (see §4.5). Figure Bl shows the cyclic magnitude and 
colour behaviour of the nine stars in the (V — /,V^)-diagram. 
In every case one observes the well known trend of bluest 
colour when the star is brightest (i.e., when phase equals 
0.0). The largest colour range is for NQ Lyr and V894 Cyg 
(the faintest ASAS star in the sample), and the smallest 
colour range is for the brightest star in the sample, NR Lyr. 
The apparent red mean colour for NR Lyr probably is due 
to its relatively large reddening (see column 5 of Table 5). 
The 'wiggles' are artificial, a result of noise and imperfect fit- 
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Figure Al. (Top) Light curve for KIC 7021124, phased with the 
period 0.6224925(7) d and tO = 54953.0 (the mean Kp magnitude, 
13.550, has not been added to the ordinate). (Bottom) Fourier 
amplitude spectrum of the detrended Ql long-cadence data, after 
prewhitening with the above period. The pattern seen here is 
similar to that seen for V350 Lyr in fig. 7 of Benko et al. (2010). 
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Figure A2. KIC 7021124 (and V350 Lyr) in the Petersen (P2/P0 
vs. Pq) diagram. The mass inferred from this diagram is ~0.75 
A4q, higher than both Al(puls) and A4(evoI). 
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Figure Bl. Looping behaviour in the H-R diagram for the nine Kepler non-Blazhko RR Lyr stars with V, I photometric data in the 
ASAS-North catalog. The stars are ordered according to increasing period, and the points along the mean light curves occur at every 
1/lOOth of the phase. Total risetimes, V amplitudes, {V) magnitudes and {V — I) colours for each star are given in the last four columns 
of Table 2. 



ting of the mean light curves; nevertheless, the mean trends, 
magnitudes and colours should be significant. 



